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ABSTRACT 


Some new developments relevant to the design of single-element 
airfoils using potential flow methods are presented. In particular, 
the role played by the non-dimensional trailing edge velocity in 
design is considered and the relationship between the specified 
value and the resulting airfoil geometry is explored. In addition, 
the ramifications of the unbounded trailing edge pressure gradients 
generally present in the potential flow solution of the flow over an 
airfoil are examined, and the conditions necessary to obtain a class 
of airfoils having finite trailing edge pressure gradients developed. 
The Incorporation of these conditions into the inverse method of 
Eppler is presented and the modified scheme employed to generate a 
number of airfoils for consideration. The detailed viscous analysis 
of airfoils having finite trailing edge pressure gradients demonstrates 
a reduction in the strong Invlscld-vlscld interactions generally 
present near the trailing edge of an airfoil. 
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CHAPTER I 


INTRODUCTION 


Low-Speed, Single-Eleroent Airfoil Design 

Much of the current research effort applied to low-speed airfoils is 
directed toward the analysis and design of raulti-eleraent sections, which 
incorporate high-lift devices such as multiple-slotted flaps and movable 
leading edge slats. The use of such airfoils permits a broad range of 
performance through the integration of an airfoil that is suitable for high- 
speed cruise with a configuration capable of high-lift for take-off and 
landing. In spite of the strong interest in the complicated flow phenomena 
connected with the multi-element design, considerable motivation remains for 
the study of single-clement wing sections. For example, a number of practical 
applications exist, including low-speed recreational aircraft, sailplanes, 
helicopter and windmill rotors, and aircraft in the expanding arena of 
remotely piloted vehicles (RPV's), for which either the cost and mechanical 
complexity of a high-lift system cannot be justified, or such a section is 
unnecessary in that a broad performance envelope is not a design 
requirement. Furthermore, much of the increased understanding resulting from 
the study of single-element airfoils is directly applicable to the individual 
components of the multi-element designs. 

As detailed by various authors, most notably Wortmann [l]-[4], Eppler 
[3]-(7], and Miley [8], the modern methodology of low-speed airfoil design 
involves relating the aerodynamic performance features sought to particular 
characteristics of the boundary layer and, lu turn, specifying the velocity 
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(pressure) distribution around the airfoil necessary to achieve those boundary 
layer characteristics which give rise to the desired performance. Once the 
required velocity distribution has been established, the airfoil is obtained 
by means of any one of a number of inverse (design) procedures such as the 
exact potential flow methods, based on complex functions, developed by 
Lighthill (9), Eppler 15) and Arlinger [lUj. Whereas in the direct (analysis) 
problem every airfoil has a corresponding velocity distribution, the inverse 
problem is complicated by the fact that every velocity distribution does not 
necessarily have a corresponding airfoil. Thus, in the design process it is 
necessary to allow some flexibility in the prescribed velocity distirbution 
such that a physically reasonable airfoil can be obtained. Also, in some 
Inverse methods, including those of Lighthill (9J and Eppler (5J, the velocity 
distribution is specified through parameters which are rather indirect. 
Consequently, these methods usually require some amount of iteration in order 
to determine what parameter values actually achieve the desired velocity 
distribution. 

One particular low-speed airfoil research subject which has received 
considerable attention is concerned with the theoretically interesting 
question of how much lift can be generated by a single-element airfoil without 
using active means of boundary layer control. In order to make the problem 
more tractable, most efforts In the area have considered only flows which are 
fully attached. Typically, airfoils in application achieve maximum lift when 
the lift increase due to an increase in angle of attack is Just equalled by 
the lift losses due to separation; and consequently, the production of high 
lift Is usually accompanied by relatively high drag* Thus, by limiting the 
maximum lift airfoil problem to fully attached flows, not only are the 



compl LcaCions ot analyzing separated flows eliminated, but the desirable 
result that the airfoils generated attain high lift with comparatively low 
drag is also achieved. 
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As discussed by Smith [llj and Llebeck [12], the velocity distribution 
formulated for the purposes of maximizing the lift generated by a single- 
element airfoil having fully attached flow is ot the form shown in Figure 1. 

On the lower surface, the desired velocity distribution is simply that which 
is as close to stagnation over as much of the lower surface as is possible. 

On the upper surface, it is dictated that the flow accelerate rapidly from the 
leading edge stagnation point to a level of constant velocity (rooftop 
velocity) that is to be maintained as far aft as possible and still permit 
pressure recovery over the rear of the airfoil without introducing flow 
separations due to an excessive adverse pressure gradient. As first adopted 
for airfoils by Liebeck and Orrasbee (13), the theoretical recovery that allows 
the longest rooftop by achieving the recovery of a given pressure in the 
shortest distances possible is the zero skin friction approximation of 
Stratford [lA]. The notion of zero skin friction implies that separation is 
everywhere imminent along the Stratford distribution. The ideal maximum lift 
velocity distribution was further defined by Orrasbee and Chen (15) In that the 
optimum relationship between the maximum velocity on the upper surface, the 
rooftop velocity, and the trailing edge velocity was determined. This 
relationship, however, does not specify the magnitude of either the rooftop 
velocity or the trailing edge velocity, but only the optimum ratio between the 
two* As a consequence, disregarding the small adjustment in rooftop length 
necessary because the local Reynolds number in that region changes, increasing 
the value of the trailing edge to free-stream velocity ratio, V^^/U, allows 
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the level of the entire upper surface velocity distribution to increase while 
the margin against separation offered by the Stratford distribution is 
unchanged. Thus, increasing V^g/U Is an extremely effective means of 
increasing the amount of lift generated by a particular design. At this 
point, however, although there is no doubt that for practical airfoils the 
value of Vj£/U must have an upper bound, it is not evident what that bound is. 
Because of its strong influence on the amount of lift generated, and because 

its specification can be connected to minimizing drag through the attainment 

of fully attached flow, it would be of some benefit to low-speed airfoil 
design to better understand the influences that the specification of Vj^/U has 
on the aerodynamic characteristics of an airfoil, as well as how the 
specification physically Impacts the geometry of the profile which results. 

Consideration of Flow Conditions in the Vicinity of the Trailing Edge 

In surveying the available literature to better understand what limits 

exist on the value of V^g/U, it is found that conventional airfoils, such as 

those catalogued by Abbott and von Doenhoff [16], exhibit a trailing edge 
velocity in the neighborhood of eight to nine-tenths of the free-stream 
value. Liebeck [12], indicates that the bounding value for Vjg/U in potential 
flow is unity. This conclusion is advanced through the argument that for a 
ciisped trailing edge, symmetrical airfoil at zero angle of attack, the 
trailing edge velocity approaches that of the free-stream as the thickness of 
the airfoil goes to zero, i.e., a flat plate. Thus, it is supposed that any 
thickness or lift generation requires that the trailing edge velocity bt Isss 
than free-stream. While no disagreement is taken with these notions, which 
can be demonstrated using the familiar Joukowsky transformation. It remains to 



develop them in a more general manner and to indicate how the phy5?ical 
characteristics of an airfoil are affected by the specified value of 
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Another characteristic of the flow in the vicinity of the trailing edge 
which warrants further consideration, and will be developed in detail later, 
is the fact that, in general, the complex pressure gradient, p^^-ipy, at the 
trailing edge is infinite. That p^^^ is infinite is readily confirmed by 
noting the vertical slope that occurs at the trailing edge in the pressure 
distributions of Joukowsky airfoils, for example. Although not as readily 
observable, the value of p^^ exhibits similar behavior. 

In addition to the presence of infinite trailing edge pressure gradients, 
there are numerous examples in the literature, including References [17]-[19], 
of airfoils In which the potential flow velocity distribution is characterized 
by a large velocity differential between the upper and lower surfaces over the 
aft portions of the airfoil and in the region of the trailing edge* As is 
demonstrated by Figure 2, this large velocity differential Is introduced as it 
increases the area enclosed by the velocity distribution and the resulting 
Increase In aft pressure loading on the airfoil manifests itself through 
increased lift production. While It might appear, in some cases* that the 
velocity distribution is closed by a vertical slope connecting the upper and 
lower surface velocities, because potential flow theory requires that each 
point In Che flow field have a unique velocity then, as discussed by Nonweiler 
(20], there can be no difference between the upper and lower surface flow 
velocities at the point where they meet and flow into the wake* Thus, 
although Che recovery distributions of a number of airfoils appearing in the 
literature, including Chat of Figure 2, have been formulated using a value of 
Vt^/U in excess of unity, in actuality, the value used is not that of the 
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trailing edge hut, rather, corresponds to a portion on the upper surface 
slightly upstream of the trailing edge. From this point, the fluid Is 
decelerated very rapidly to the actual velocity at the trailing edge. This 
value is in common with chat resulting from accelerating the flow on the lower 
surface in the immediate vicinity of Che trailing edge through a steep, 
favorable pressure gradient. While it is quite clear that the viscous effects 
prevent the full realization of the lift gains predicted by potential flow 
methods, exper results have indicated, in some applications, that the 

proper implementation of large velocity differentials between the upper and 
lower surfaces to very near the trailing edge can be of some benefit. 

An interesting example of an airfoil designed to exploit the benefits of 
a large amount of aft loading is that of Kennedy and Marsden [18J. The 
potential flow analysis of this airfoil, shown in Figure 2, yields a lift 
coefficient of 3.81, resulting from the design velocity distribution based on 
a upper surface trailing edge velocity of 1.2U. While experimentally, the 
lift coefficient at the design conditions was found to be an impressive 2.64, 
it was obtained at the expense of relatively high drag resulting in a maximum 
lift-to*drag ratio notably less than those obtained for other high lift 
designs* It was also found experimentally that viscous effects reduced the 
upper surface velocity just upstream of the trailing edge from its potential 
flow design value to 1.07U. To further demonstrate the viscous Influences, 
Figure 3 compares an off-deslgn potential flow velocity distribution of the 
Kennedy and Marsden airfoil with one obtained experimentally, Reference 
[18]« In addition, the figure includes results obtained using the CKUMFOIL 
code developed by Mead and Melnik [21]* In the classical method of correcting 
Invlscid flow results for the effects of viscosity, the displacement thickness 
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of the boundary layer on the airfoil is calculated and added to the orlj^lnal 
profile. This results in an equivalent body which can be analyzed to 
approximately account for the influence of the boundary layer. This 
procedure, however, ignores additional viscid-lnvlscid interactions, each of 
which has an influence on the inviscld result equal to that of the 
displacement thickness. The GRUMFOIL code remedies the diflciencies of 
previous methods for the analyses of the flow over an airfoil by Incorporating 
a complete Interacting boundary layer formulation which, in consideration of 
application to low-speed flows. Includes the effects of the boundary layer 
displacement thickness on the airfoil and In the wake, wake curvature effects 
arising from the turning of low momentum fluid In the wake along curved 
streamlines, and the effects of strong viscous Interactions in the vicinity of 
the trailing edge. As Figure 3 demonstrates, the potential flow velocity 
distribution over the airfoil is modified considerably by viscous effects. In 
particular, the steep trailing edge velocity gradient is rounded off to the 
extent that the actual trailing edge velocity of .94U is easily Identified. 

It is also evident from Figure 3 that, because the viscous Influences are 
amplified In high lift situations, the conscientious design of high lift 
airfoils for actual applications must involve the integration of viscous flow 
analysis with the potential flow design method. 

In the potential flow inverse scheme of Eppler (5), an attempt is made to 
achieve the lift-increasing benefits of a large upper and lower surface 
velocity differential near the trailing edge and also, in a somewhat 
qualitative manner, account for Che viscous effects. While the potential flow 
formulation requires that Che upper and lower surfaces have equal velocities 
at the trailing edge, the manner in which geometrical closure of the airfoil 
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is achieved allows for large velocity differences between the upper and lower 
surfaces to a point very near the trailing edge* Consequently, this trailing 
edge treatment allows for velocity distributions which, in the /iclnity of the 
trailing edge, are very similar to those which occur when viscous influences 
are considered. While this model is :learly an improvement on the usage of 
the nearly infinite trailing edge gradie* . the question arises as to what is 
the permissible extent of recovery before the steepened portion of the 
distribution causes flow separation problems. For example, if a Stratford 
distribution is specified for recovery on the upper surface, then, even in the 
best of situations, separation would occur where the velocity gradient 
steepens due to the closure contribution. As is demonstrated by both the 
experimental and viscous analysis results shown in Figure 3, the loss of 
loading near the trailing edge that occurs when viscous effects are considered 
generally causes the overall recovery gradients to be somewhat steeper thnn 
indicated by inviscid results. Thus, it is a possibility that the steeper 
gradients will cause severe separation problems. Even in the case of more 
practical airfoils, not pushing recovery limits as does a Stratford 
distribution, only the gentlest distributions will have sufficient momentum in 
the boundary layer tc overcome the very steep adverse pressure gradient 
Introduced in the vicinity of the trailing edge by the presence of large 
amounts of aft loading* To reduce these problems, Eppler (S] discusses the 
fact that, while it is desirable to reduce the upper surface adverse pressure 
gradient in the closure region as much as possible by increasing the value of 
the trailing edge velocity ratio, such relief is limiced* As the value of 
Vj^/U is Increased, a point is reached for which any additional Increase 
causes the upper and lower surfaces of the alrfcll Lo intersect one another 
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ahead of Che trailing edge* In addition, Nonweiler [20] notes that even if it 
were possible to completely eliminate the upper surface adverse pressure 
gradient due to closure by accelerating Che lower surface flow through a 
strong favorable gradient and discharging the flows from both surfaces into 
the wake at a velocity corresponding to a smooth continuation of the upper 
surface recovery distribution that, as this dumping velocity would be greater 
than ch.^t of the free stream, some deceleration of the flow in the wake would 
be required. This situation is also likely to have an upstream influence 
manifested as separation on the upper surface ahead of the trailing edge. 

A demonstration of the problems that might occur using the trailing edge 
treatment of Eppler [5] is provided by the high-lift, single-design-point 
airfoils of Thompson [I9j. These airfoils make use of the Stratfo"-d recovery 
distribution on the upper surface. As shown in Figure 4, the closure 
contribution at the trailing ed^e of these airfoils extends over such a narrow 
portion of chord length that, like the Kennedy and Marsden airfoil of Figure 
2, it appears as though the large velocity difference between the upper and 
lower surfaces is adjusted at the trailing edge through an infinite velocity 
gradient. As illustrated in the insets of Figure 4, however, it is found that 
the value of V^g/U used to generate the upper surface recovery distribution is 
actually an upstream point at which is initiated a very rapid deceleration of 
the flow around the small protuberance located at the rear of the airfoil. 
Thus, the actual Vjg/U is considerably less than the value in excess of unity 
which was employed in formulating the Stratford distribution. Consequently, 
rather than the desired velocity distribution, which includes recovery through 
a Stratford distribution extending to the wake, the airfoil generated has 
recovery through a Stratford distribution followed by a violent recorapression 



at a point near the trailing edge. Not only might this result lead to flow 
separation problems and a thicker wake than is necessary, but more 
importantly, the discrepancies that exist between the intended velocity 
distribution and that obtained introduces a degree of uncertainty into the 
design procedure which is undesirable. While a part of these problems can be 
attributed to erroneous interpretations the human-computer iteration 
process which is necessary with the inverse method used for these designs, a 
more significant element was the fact that the values of Vj^/U used in 
formulating the desired velocity distributions were much too high- While this 
possibility was clearly acknowledged, it was also emphasized by Thompson [19] 
that Che literature contained very little information to aid in the selection 
of this Important parameter. 

Objectives of the Present Investigation 

As part of the desired velocity distribution, the value of the trailing 
edge velocity ratio must be specified in most airfoil design procedures; 
however, there is little information in the existing literature Co guide in 
its choice. Consequently, as has been discussed by several authors, including 
McMasters and Henderson [22], V^g/U is perhaps one of the most difficult of 
the required parameters to determine. While in many design exercises, such as 
those directed at low drag, the situation is less critical in that the desired 
aerodynamic characteristics are less sensitive to an optimum selection of 
V»j.g/U, the difficulty is heightened in the maximum lift problem because of the 
strong dependency of the lift generated to the value specified. Thus, in 
order to provide some guidance for selecting the value of design 

of low-speed, high-lift airfoils, the first part of the present investigation 
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considers the relationship between the maximum values of obtainable and 

profile characteristics such as thickness and camber. 

Considering the modification of the potential flow maximum lift velocity 
distributions by viscous effects, there is clearly a trade-off between the 
lift gained by maintaining a large velocity differential between the upper and 
Lower surfaces in the vicinity of the trailing edge and the lift lost by 
separation. While there are examples of potential flow designed airfoils that 
achieve some portion of their design goals, Kennedy and Marsden (18] and 
Sivier, et. ai • (23), there are others, such as the designs of Thompson (19] 
experimentally Investigated by Moore (24), which exhibit extremely poor 
performance attributed primarily tc widespread flow separation. In the 
context of designing airfoils having predictable characteristics, the 
inconsistency of results for airfoils which make use of large velocity 
differentials near the trailing edge leaves much to be desired. Thus, until 
additional guidance is available, the design process is limited in not knowing 
precisely to what extent the potential flow results in the vie nity of the 
trailing edge will be modified by viscous interactions. In order to eliminate 
some of the uncertainty in using potential flow airfoil design methods, the 
second part of this research is directed toward the development of a class of 
airfoils in which the viscous interactions in the trailing edge region are 
minimized by the introduction of a condition to insure that the pressure 
gradients at the trailing edge are bounded* Not only does the enforcement of 
this condition permit the fluid on the airfoil to flow smoothly Into the wake; 
but in addition, it typically excludes the possibility of the closure 
contribution causing steep upper surface adverse pressure gradients in the 
vicinity of the trailing edge* The practical implementation of this condition 







The final phase of the research effort to be presented is directed toward 
exploring the influences on the geometry and aerodynamic characteristics 
caused by the Introduction of the finite trailing edge pressure gradients 
condition into the design process. To facilitate this exploration, a number 
of airfoils having finite trailing edge pressure gradients are presented and 
compared to airfoils which are as similar as possible but having unbounded 
trailing edge pressure gradients. 
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CHAPTER II 


MAXIMUM TRAILING EDGE VELOCITIES 


Preliminary Remarks 

As Che ratio of the trailing edge to free-sCream velocity, Vf£/U, has a 
global Influence on Che lifting capability of a given airfoil, it Is important 
:o determine Che maximum value that this ratio can attain and how its 
specification affects other airfoil design considerations. While, as it has 
been noted, there are assertions in the literature regarding a potential flow 
upper bound of V.p£/U equal to unity for a flat plate at zero angle of attack. 
Reference [12] for example, these are apparently based on results formulated 
through Che use of the Joukowsky transformation. A survey of Che literature 
has not revealed an extension of this conclusion for the more general 
transformation mapping a circle into an airfoil. In order Co develop such a 
demonsCrati ui, a series of truncations of the general transformation, the so- 
called von Nises transformations, will be used to show that the addition of a 
finite number of extra terms to Che series does not alter the flat plate 
results. While this does not result in a general proof concerning Che 
ultimate value chat Vf£/U can attain, it does indicate that Che upper bound 
for Che v^ry broad class of von Mises airfoils is indeed unity. More 
importantly, however, the developMnt provides some insight into how the 
specification of V^£/U influences the geometry of an airfoil. 

In the course of numerical design studies in which it was attempted to 
raise V.^£/U to as high of value as possible, it was found that while Che 
potential flow formulation of Che problem places no restrictions on the upper 
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bound of V.£.£/U, the resulting airfoils are physically restricted in that as 
increases are nade in V»pg/U, a point is eventually reached for which further 
increases cause the upper and lower surface to cross one another such that a 
profile .laving a region of negative thickness is generated. Thus, this 
situation indicates that an additional constraint must be imposed such that 
the problem becomes that of finding the largest value of for which a 

physically realizable airfoil is possible. 


Transformation of a Circle into an Airfoil 

The method of conformal mapping may be employed to analyze the flow over 
airfoils by transforming the known flow field about a circular cylinder into 
that about an airfoil as depicted in Figure 5. In this formulation, 
t * u denotes the center of the circle, r its radius, and the function 

■ u + re^^ describes the circle. The complex potential for the flow about 
the circular cylinder with circulation T can be expressed as 


F(C) 


U(c-u)e"^“ + 


Ur ia ir (C-W) 

^ 27 


(2-1) 


and the complex velocity in Che C'plane is 




( 2 - 2 ) 


The complex velocity in the airfoil plane is then given by 




(2-3) 
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where z(^) is the mapping function which takes the circle into the airfoil. 
This function will. In general, possess a certain number of critical points 
which are defined by the solution of the equation (dz/dc) * 0. A profile 
having a sharp trailing edge requires that one of these critical points be on 
the circle while all of the others are contained In the Interior of the 
circle. Thus, the critical point on the circle, s.j., maps to the trailing edge 
of the airfoil while the regions external to the circle and the airfoil are 
everywhere conformal. 

The amount of circulation present In equation (2-1) may be determined by 
introducing the Kutta condition, which requires that the flow velocity at the 
trailing edge be finite and continuous. From equation (2-3), the complex 
velocity at the trailing edge Is given 

w^ ■ w(z^) ■ w(c^)/z^ (2-4) 

Because a sharp trailing edge requires that be zero, the requirement that 

Wj be finite necessitates that w(c^) also be zero. The circulation which 

*16 

fixes this stagnation point at Cj * re 4- p is found from equation (2-2) to 
be 


r - 4xRU sin (a + S) (2-5) 

As presented in Keramcheci [26], a non-zero trailing edge closure angle 
requires a stagnation point at the trailing edge of the airfoil. Thus, In 
Investigating upper bounds of V^g/U, it Is necessary to consider airfoils 
whose trailing edges close in a cusp and thereby have non-zero trailing edge 
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velocities. In order to obtain an expression for the trailing edge velocity 
of this class of airfoils, L'Hopital*s rule is applied to the indeterminate 
form obtained from equation (2-4) such that 


|w(z^)| * lira 


'(OL 


( 2 - 6 ) 


and, using equations (2-2) and (2-5) to obtain 


, 2U e 

^ m I. .1 

T r 


cos (o + 6) 


(2-7) 


the trailing edge to free-stream velocity ratio is found to be given by 


TE _ 2 |co8(g B) 


l*x 


(2-8) 


For the Joukowsky transfomatlon. 




(2-9) 


where, as can he verified by consideration of Figure 5, 


cos 6 


(2-10) 


the trailing edge velocity ratio la then given by 




(2-U) 
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As shown in Reference [26], as well as in any one of a number of classical 
aerodynamic texts, the thickness of the Joukowsky airfoil is increased as the 
magnitude of e is increased, *«’hile increases in B increase the amount of 
camber. Thus, from equation (2^11), it is clear that for the case of the 
Joukowsky transformation, the presence of any thickness or camber reduces 
Vig/U from the value of unity which exists for the flat plate at zero angle of 
attack. 

The Von Mises Transformation 

The general transformation which maps a circle into an airfoil is 
expressed as 


z ■ z(c) 


m 

c > I 

n-l 



( 2 - 12 ) 


where, as can be seen in Figure 5, r ■ x + ly, C ■ C + In , and, in general, 
Che c„'s are complex quanclcles. Note chat Che mapping funcClon, z(c), 
aacisfies Che requiremenC of noC alCering Che flow field aC infiniCy. 

Now, consider a CruncaCed form of Che general cransformacion which, afcer 
dlfferenciacion, is mriccen as 


di 

dC 


N 


I 

n-l 


n c 

n 



(2-13) 


This Cransformacion muse noc possess any singular poinCs outside the 
generating circle, although one singular point, denoted by must be on the 

circle. Since the remaining singular points, given by r,^, are 

zeros of equation (2-13), in factored form, that equation becomes 
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dz 

dC 


(I (1 (1 [1 -/) 


(2-14) 


where Che zero locations are Indicaced by 


C - C + ln„ (2-15) 

n n n 

Ac this point, it la seen chat Che mapping function may be generated by 
postulating the locations of Che N zeros which lie, along with the origin, 
within a circle of radius r in the (-plane* Without loss of generality, the 
zero on the circle, which transforms to the trailing edge of Che airfoil, 
may be fixed at Che location (1,0). By choosing different sets of zeros and 
different generating circles, a great deal of flexibility exists in the shapes 
of airfoils which can be obtained. Profiles which are generated by these 
means are known as von Nlses airfoils and the mappings obtained by expanding 
equation (2-14) and integrating the result are kno%m as von Mlses 
transformations. 

The coefficients of *:he transformstion, the c^'St can be related to the 
zero locations, the by expanding equation (2-14) and equating that 

result with equation (2-13). In this manned, it is found Chat Che coefficient 
of the ( ^ term must vanish. Hence, the relation 

N 

Cp ♦ I C- • 0 (2-16) 

^ n-2 ** 

indicates that the origin of the (-plane must be Che centroid of the zeros 
located at (^, ( 2 » C^t ••• 


(^. The coefficients for von Mises transformations 
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of up to six terms, N < 6, are related to the zero locations as follows: 

Cj - ic/ - * ^U* ^5* ^ '^S ^ '^6^ 

^4(^5 ■•■ (2-18) 

C 2 - j {CT.lt2(C3 ♦ «4 ^ ^ ^^6^ * ‘^3^«4 ^ S ^ ^ 

4 C^(Cr > C^) + CjCj + C2t'3^'4 ^ '5 ^ 

+ 54(^5 + Cg) + ^ ^ *^5^61 ^ 

C3 - - j 1 Ct[«2<^3«4 * '3'^5 * ^ ^ 

* «3<^4<^5 "■ h ^6 * * '4«5'6l ^ '^2^'3^'^4‘^5 

^ '^4«6 ^ * «3'=4'5'6^ 

^4 " i i^T^^2^'^3^4S ■** ^3^4 ^6 ■*■ '^3^s''6 '*’ ^4^^6^ 

♦ ?3«4«5'^6l ^ Wkhh^ <2-21) 

‘^S - “ T WW4S'’6> <2-22) 
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For a two term transformation, N«2, equation (2'-17) Indicates that 
C 2 ” Cj snd equation (2‘-14) reduces to the familiar Joukowsky case, 


— ■ f 1 - —1 

dC 2 ^ 


(2-23) 


Prom equation (2-13), and the Joukowsky transformation is obtained In 

the form of equation (2-12) as 


; + 


?t2 


(2-24) 


While Joukowsky airfoils are limited to a circular arc camberline and a 
maximum thickness at approximately the quarter-chord position, such 
restrictive limitations do not exist for von Mises airfoils. By appropriately 
locating the generating zeros and the origin of the transformation circle, it 
Is possible to approximate a desired shape through an extreswly wide range of 
airfoil geomatries. While the zero locations govern the basic thickness 
distribution and camberline shape, the center location of the generating 
circle can be used to influence the overall amounts of thickness and camber as 
in the case of Joukowsky airfoils. Displacements of the circle in the 
direction of the real axis, Ug, primarily influences the thickness of von 
Hiaes airfoils while displacements in the direction of the imaginary axis, u^, 
the camber. 

A few exa^tles of the profile variations possible using the von Mises 
transformation are provided by the airfoils shown in Figures 6 and 7, 
generated uelng four term transformations, and those shown in Figures 8 and 9, 
obtained using six terms. The sero locations used to generate these airfoils 
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and Che resulting cransfornacion ccetficiencs are Indicated In the figures 
with the complex number pairs denoted by the Z^'s and Che 
respectively. The real and imaginary components of the origin of the 
generating circle are given in the figures, respectively, as UR and Ul. 

Because the von Hises transformations provide for such a wide range of 
possibilities, and airfoils of practical interest can generally be 
characterized by relatively few inflections in Che thickness distribution and 
camberline shape, the conclusions developed using the von Mises 
transformations are considered to be a,>plicable over a broad range of 
airfoils. Tr, addition, Che development is conducted using a relatively small 
number of terms in the transformation because, as demonstrated by the examples 
presented, large variatons in profile geoMtries are possible without a large 
number of terms. This is further Justified by the fact that the contributions 
made to the resulting airfoil shape by terms ot increasing order rapidly 
become insignificant compared to the Influence of the first few terms. 

Maximum Trailing Edge Velocity Ratios for Physically Realisable Airfoils 

From equation (2~f), it is seen that V^£/U is inversely proportional 
|sj’*| and r. Thus, the maximum value of V^/U occurs when Isj*' Is minimized 
and r is fixad at a minimum valua for the given Iz-f”! which insures that the 
airfoil generatad is free of any regions of negative thickness* Denoting the 
complex coefficients of the transforsMtion function as ■ a„ t lb,,, taking 
Cj ■ O*0)» and performing the necessary operations on the truncated form of 
equation (2-12), |sf“| is found to be 



N-» , N-l .. 

IV'I • ^ n(n+l) ( I n(n-M )bj^] 

n»l n"l 


( 2 - 25 ) 


Because Che second tern or Che righc of Chls expression Is poslcive, Che 
mlnlaum value of |z<f‘*|, for a given seC of a^'s, occurs when Che second Cerm 
zero. While chere are non-zero values of Che chac can achieve Chls, 

Che second Cerm always disappears for synmeCrlcal airfoils which are generaced 
when and all of Che bj^'s are zero. Consequencly, as developed lurcher In 
Appendix A, for deCermining Che maximum value of V^^/U pos. >le for an airfoil 
generaced by a fixed sec of Ic Is sufficienc co consider only Che 

symmecrlcal profile which csn resulc. While chere mey be non-symmecrical 
secclons having a value of V^£/U as large, chere csn be none for which Che 
value is larger. By eiiminsCing che need Co consider oCher chan symowctical 
airfoils, Che decerminscion of che Mximum accainable values of for 

physically realizable airfoils of che von Miscs family is significancly 
simplified in ChsC Che chickness disCribuCions, which muse exclude negacivc 
ehicknesses, art dependenc only upon and Che a^'s. Through Chase 
considerations, che problem of deCermining che conditions insuring chac a 
cransformacion yield an airfoil which does not cross-over itself is reduced to 
chac of finding che minimum value of r for which che upper surface coordinates 
of Che generaced profile are all non-negacive. 

To proceed in chls development, che representation of the generating 
circle in che complex C'Tlane, as previously noted, is given by 


C • C ♦ in ■ re^^ ♦ w 


( 2 - 26 ) 
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Thus, the coordinates of the circle are 

C = rcos ^ 
n * rsin4> + 

where, with * (lfO)» it is evident from Figure 5 that 

} 1 I r / 1 , 2 *1 1/2 

r « 1 1 - M I * ^ i 

Setting to zero and defining 

f » (I - COS(>) 

allows the preceding relationships to be written as 

r - 1 - Ur 

C - 1 - rf 
n • rsln^ 


(2-27) 

(2-28) 


(2-29) 


(2-30) 


(2-31) 


(2-32) 


(2-33) 


The relationship between the real part of the transformation 
coefficients, the a^'s, and the radius of the generating circle Insuring that 
the symmetrical airfoil developed Is characterized by positive thickness 
everywhere has been analytically determined for von Mlses transformations of 
up to six terms* Before this case Is considered, however, it is instructive 
to examine Che derivation with only three terms as the key points of Che 
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development are retained while the algebra is simplified considerably. 

The three-term von Mises transformation having real coefficients is 
written as 


a 


(2-34) 


For the symmetrical airfoil generated to have non-negative thickness 
everywhere, it is required that 


Im(z) ■ y > 0 


(2-35) 


for 0 < ^ < V. Substitution of equation (2-26) into (2-34), and making use of 
the condition given by equation (2-33), yields 


a^n 2^2^^ 


(2-36) 


which may be rearranged to give 


+ n^) - 2 b 2 ^ > 0 


(2-37) 


Using equations (2-32) and (2-33), the quantities required in equation (2-37) 
are found to be: 


5^ . p2^2 _ 2rf + 1 


(2-33) 


2 2 2 
■ r sln^t 


(2-39) 
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+ n^) » 2 rf(r-l) + 1 


(2-40) 


* 4rf (r-l)[rf (r-1) + l] + 1 


(2-41) 


Substitution of the above quantities into equation (2-37) gives 


4rf(r-l)[rf(r-l) + l] + 1 - aj[2rf(r-l) + 


- 2a, [l - rf] > 0 


(2-42) 


The terms in this expression which are independent of r and f, that is, 1-a^- 
2 a 2 « may be eliminated by making use of the Kutta condition which requires, 
for ttie general case, chat 


N-1 

*t' • ‘ ^ " *n • ° 

n*l 


(2-43) 


For N~3, this expression indicates that 


1 - a| - 2a2 ■ 0 


(2-44) 


2r(r-l)*f + 2(r-l) - aj(r-l) + a 2 > 0 


(t) 


Applying this result and factoring rf from what remains, equation (2-42) 
becomes 


P 


(2-45) 
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As the first term in this expression is non-negative, given that r is non- 
negative, the most critical situation for meeting the condition of the 
inequality occurs when f ■ 0, i.e., when cos<^ » 1. For this case, equation 
(2-45) becomes 



4(r-I) - 2a^(r-l) + 2a2 > 0 


(2-46) 


Using this result, it is found that in order to prevent the airfoil obtained 
through the mapping function from having regions of negative thickness, it is 
necessary to require that 


r > 



(2-47) 


In addition, the value of equation (2-25), is given by the 

expression 


I*t"I “ I + eaj I 


(2-48) 


Using equation (2-8), the maximum value of ^'te /U which can be obtained from a 
given three-term transformation is 



2 COSO 



2(2-aj)cosa 


(2-49) 


For this case, the maximum V.j.g/U clearly occurs when the angle of attack, a. 
Is zero* To determine the transformation which affords the largest V<p£/U 
overall, equations (2-18) and (2-19) are used to relate the coefficients to 
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the locations of the generating zeros. Thus, for symmetrical airfoils 
generated with all of the generating zeros located on the real axis such tha' 
the are all zero, the coefficients are given by 

a^ - 1 - (2-50) 

Substitution of these relations into equations (2-49) and setting a equal to 
zero yields 



(1 V3)' 


and using equation (2-17), it is found that 


C3 - - (1 + 


Equation (2-32) thus becomes 



4(C2^ + + 1) 

(^ 2 ^ + ?2 - 2 )^ 


(2-52) 


(2-53) 


(2-54) 


Maximizing V^g/U with respect to C2 indicates that the overall largest value 
of V^g/U occurs when ?2 ■ "1 and C3 ■ 0 or * 0 and Cj • - 1 . While it Is 
also possible to generate symmetrical airfoils with non-zero values of the 
n^'s, for example, when C2 and C3 * ~*5 and ^2 " ~ ^3, it is found that the 
value of V<fg/U is maximum when h2 * **13 ■ 0 . Thus, for determining the 
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bounding values of V^g/U, it is sufficient to consider the case of having 
the n^*s set to zero. Thus, for the three-term von Mises transformation, the 
maximum occurs when the third terra is zero and the remaining terms are 

equivalent to those of the Joukowsky transformation which yield the flat plate 
airfoil. 

Now consider the case of the six-term von Mises transformations having 
real coefficients, as given by 


z 





a 

n 



(2-55) 


A physically realizable symmetrical airfoil is insured provided that 

5 a a - in)" 

Im{C + in + I -A; 5 } > 0 (2-56) 

n-l + n^)" 

Proceeding as in the case of the three term transformation, this equation is 
expanded to obtain 



^ 2>,5 
+ n J 


» n^)‘ - * .3(n^ - 


* 4«4C(n^ - E^)(«^ + n^l ♦ ajClOE^n^ - 55“ - n‘) > 0 (2-57) 


After making the proper substitutions and expanding, the terms in the 
resulting equation which are independent of r and f are again eliminated by 
use of the Kutta condition, equation (2-43). Factoring 2rf from the remaining 
expression yields an inequality of the form 
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4 3 2 

Af + cr 


+ Of + E > 0 


(2-58) 


where expressions for the coefficients are given in Appendix B. For 
0 < < TT, f, given by equation (2-30), is never negative. Hence, when A, B, 

C, and D are non-negative, the inequality will be satisfied as f goes to zero 
provided that E is non-negative. Additionally, for the case when A, B, C, and 
0 are not all non-negative, it has been demonstrated numerically that 
requiring E to be non-negative is still sufficient to guarantee that the 
condition of the inequality is net. From these results, the attainment of a 
non-negative thickness distribution on an airfoil resulting from a six-term 
von Mises transformation requires that 


5 - 4aj - 7^2 “ ^^3 * ^^^4 * ^0^5 
r > 5 - 4.J - 6.2 - 5.3 + lOa^ 


(2-59) 


It should be noted that the expression degenerates into the appropriate forms 
for transformations having fewer than six terms. 

The results of maximizing for transformations of up to six terms in 

a manner analogous to that used for three terms suggests that for * (1,0), 
a specified ^2 location, and any number of remaining zeros, say k, that 
achieves a maximum when the k zeros are all positioned at the location defined 
by 




0 


(2-60) 


Furthermore, Che value of V^g/U la found to Increase as Che value of ^2 
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approaches -1.0 from either direction. When ^2 equal to -l.O, the maximum 
value of V.pg/U reaches unity and requires that ail of the other zeros be 
located at » 0 and consequently, have no Influence. Again, this indicates 
that the largest possible value of V.j.g/U for a physically reasonable airfoil 
is unity and occurs when the profile is a flat plate at zero angle of attack. 

Numerical Results and Discussion 

In order to simplify the determination of the minimum generating circle 
radius for which an allowable airfoil is obtained, the development in the 
preceding section considered only syaunetrical profiles in order to uncouple 
the thickness distribution from the b^'s. Note, however, that if the airfoils 
under consideration have neither excessive camber or thickness, the resulting 
condition provides an approximation which is still useful to insure the 
generation of physically realizable airfoils. For non-symmetrical profiles, 
due to the fact that a chordwise location from an upper portion of the 
generating circle, x(9), does not correspond to the same location transformed 
from the lower portion, x(-9), the thickness distribution becomes coupled to 
the b^'s. If the airfoil is not excessively cambered, however, then x(9) will 
be approximately equal to x(-9) and consequently, the thickness distribution 
is only weakly Influenced by the non-zero b^*s. Thus, the positive thickness 
condition remains approximately correct. Similarly, in practice, it has been 
found that the condition is reasonably valid for small, non**zero values of 
and only excessively cambered profiles actually might have regions of negative 
thickness even though the positive thickness condition is satisfied. Thus, 
for most practical purposes, the condition developed is useful in generating 
physically realizable von Mises airfoils whether or not they are symmecrtcal. 
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To further understand the utility of the positive thickness condition, a 
digital computer program was written to calculate the transformation 
coefficients and airfoil coordinates resulting from an input generating circle 
and a set of zeros. As previously Indicated, on the figures to be discussed, 
the zero locations and the transformation coefficients are each listed as a 
complex number pair denoted by the and respectively. Considering 

the condition on the minimuia allowable radius such that » 1 - the 

value of URMAX presented in the figures is the maximum allowable real axis 
location of the generating circle center which results in a reasonable 
airfoil. Depending on which situation is most critical, this value is based 
on either the positive thickness condition, as determined from equation 
(2-59), or, on the requirement that all of the zeros lie within the generating 
circle. The radius of the mapping circle used to generate the airfoil shown 
is denoted in the figures as R, while the real and imaginary components of its 
center are UR and UI, respectively. The value of V^g/U for the airfoil is 
identified by VTE. 

To demonstrate the relationship of camber and thickness to the value of 

Vjg/U, a comparison of the airfoil shown in Figure 10 to that shown in Figure 

7 is informative. The symmetrical airfoil of Figure 10 has been generated by 

locating the zeros such that the are all zero, while the have 

essentially the same values as those of the section shown in Figure 7. As the 

X N-l 

transformation for the cambered airfoil yas determined so that n(n>l) b^ 
is zero, the value of Iz^'l both airfoils as can be seen from 

equation (2-*2S). Thus, with both airfoils having the same r, the values of 
Vjg/U from equation (2-8) are equal. In essence then, the airfoil of Figure 7 
has traded an appropriate amount of thickness for camber such that there is no 
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reducClon in V<|.g/U. As is also demonstrated by the expression for the value 
of Vjj;/U in the case of Joukowsky airfoils, equation (2-11), the trading of 
some amount of thickness for camber Is possible in all cases provided that 
some thickness is available to trade< Hence, in the case of the flat plate it 
is not possible to obtain a cambered plate while remains equal to 

unity. 

The effect of varying on the airfoil resulting from a given 
transformation is exemplified by Figures 9, 11 and 12. The airfoil of Figure 
9 is developed using the maximum value of as determined from the condition 
for positive thickness, which is approximately correct for the cambered 
profile shown. Thus, it can be considered that the point where the upper and 
lower surfaces cross Is located at the trailing edge of this profile. As the 
value of Is made greater than the point of cross-over moves from the 

trailing edge toward the leading edge. The airfoil shown in Figure 11, using 
the same set of generating zeros as the profile of Figure 9, exhibits this 
behavior. As the value of is increased above the value used in Figure 11, 
the point of cross-over on the resulting airfoil moves further and further 
forward. In some cases, this continues until the cross-over point reaches the 
leading edge, after which, further Increases in cause the thickness 
distribution of the resulting airfoil to be totally negative. In this 
situation, the upper parts of the generating circle are mapped to the lower 
surface of the airfoil, and the lower parts of the circle to the upper surface 
of the airfoil. An alternative possibility exists as, in some cases, before 
the cross-over point reaches the leading edge, a value of Uj^ Is used tfhich 
causes one of the generating zeros, in addition to that at the trailing edge, 
to lie on the circle. This results in an additional cusp appearing on the 
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boundary of Che airfoil* A further increase in Chen results in a zero 
becoming situated outside of the circle and the mapping of Che flow field is 
no longer conformal. 

While moving the origin of the generating circle in Che positive 
direction causes Che resulting airfoil to become thinner. Figure 12 shows Chat 
moving It in the opposite direction causes it to thicken. To the extent that 
Che positive thickness condition is an approximation for non-symnetrical 
airfoils, Che examples of figure 9, 11 and 12 demonstrate Chat generating a 
profile with equal Co results in Che thinnest uncrossed airfoil 

possible for a given sec of zeros. 

The next group of figures is included Co show quantitatively an 
additional non'*zero term in Che transformation influences the shape of the 
profile and the value of V^^/U. Figure 13 depicts Che Joukowsky flat plate 
airfoil at zero angle of attack. As Indicated in the figure, this airfoil has 
the limiting value for of unity. By moving a third generating zero a 

small distance away from the origin in either direction. Figures 14 and 15, it 
is seen that the profile thickness increases slightly while Vjg/U becomes less 
than unity. 

Figure 16 depicts an airfoil having the highest attainable value of 
using a six'term transformation and ^2 * ** specified location of the 

second generating zero* The locations of the other zeros for this airfoil are 
prescribed by equation (2-60). As before. Figure V indicates that a decrease 
in V^£/U results when the zeros are moved slightly from the maximising 
locations used to obtain the airfoil of Figure 16. 

Figures 18 and 19 present plots of the zero-lift, maximum values of Vf£/U 
and the symmetrical airfoil thickness ratios as they depend on the zero 



locaCion, ^ 2 * Figure 18 is produced using a Chree~term von Mlses 
transformation, while Figure 19 a six-term transformation. For both cases, 
the zeros other than the trailing edge zero, C,p, and the given zero, ^21 
positioned using equation (2-60) In order to maximize the value of V.p£/U. The 
symmetrical behavior about ^2 ~ ~ of the plots shown In Figure 18 Is due to 
the fact that, with only three terms, equation (2-60) requires that ^3 “ ~ 

(^2 '^ !)• Thus, the Joukowsky flat plate results that occurs when ^2 ~ ^nd 
• 0 also occurs when ^2 * 0 and C 3 “ -1. It Is Interesting to note that 
this behavior disappears, as shown in Figure 19, when additional terms are 
eoiployed In the transformation. In comparing the two figures in the region of 
^2 * If la apparent that the presence of more termr allows higher values of 
Vt£/U CO occur for a given thickness ratio. 

As an aid in airfoil design studies. Figure 20 sumsarizes the manner in 
which Che maximum, zero-lift trailing edge velocity ratios are limited by Che 
airfoil thickness ratio. It is important to note, however, chat the 
relationship shown should not be regarded as absolute potential flow limits 
as, by moving the zeros off Che real-axis, it is possible to obtain 
symmetrical profiles of larger thickness ratios which prxluce higher values of 
Vt^/U chan Chose indicated. As the profiles generating these results have 
excessively blunted or concave noses, and otherwise radical shapes, these 
cases were excluded from the revolts presented in Figure 20. 

Because Che specification of Vj^/U is necessary in a number of airfoil 
design techniques, it is anticipated that the preceding results will be 
useful. For example, in the fully laadnar airfoil designs of Sappuppo and 
Archer (27], the value of Mi^cted *'as high as possible so as 

to obtain high lift". Aa consistent with Figure 2U, the thickness ratio of 



Che symmetrical airfoil resulting from the surface singularity inverse 
procedure employed was 9.7%. On the introduction of camber, however, 
difficulties were reported in retaining a reasonable thickness ratio and that 
of the resulting airfoil was only 4.2%. While it was concluded that the 
thinness of these profiles was due to Che low Reynolds number specif icatluu 
and employment of the fully laminar Stratford recovery distribution, the 
present investigation indicates that it is, more likely, a result of the high 
value of Vjg/U specified* Furthermore, in light of the discussion regarding 
the trade-off of thickness and camber necessary to maintain a given value of 
Vjg/U, the difficulties encountered when camber was introduced are to be 
expected* 

As already noted, another area in which insight into the factors 
influencing the selection of V^/U would be most valuable is in the high lift 
airfoil design procedure detailed by Thompson (19]* In this method, the 
formulation of the Stratford recovery distribution for the uppev surface is 
initiated by seK ;.tlng the value of V^^g/U* As has been discussed, the 
literature provided little guidance for choosing reasonable values of V^/U, 
and those used by Thompson were unreasonably high* Thus, the present work 
provides the meens by which the specification of Vjg/U can be made with 
greater understanding* 
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CHAPTBR III 


AIRFOILS WITH FINITE TRAILING EDGE PRESSURE GRADIENTS 


Conditions Required for Finite Trailing Edge Preasure Gradients 

In considering the results of Interacting boundary layer theory to model 
the real flow around an airfoil, It is found that the diaplacement thickness 
of the boundary layer and of the wake increase rapidly in the vicinity of the 
trailing edge, and a discontinuity generally exists where they Join at the 
trailing edge* This behavior has been considered in numerous references, such 
as the classical papers of Preston, et. al., [29|-[31]. as well as in more 
modem sources such as Reference (32], and, as shown by Melnik, et* al*, |33j 
Is attributable to singularities which occur at the trailing edge of the 
invlscld solution* These unbounded quantities, present In the case of any 
airfoil carrying a non-zero load at the trailing edge, lead to a breakdown of 
the boundary layer approxlmatlops such that. In the vicinity of the trailing 
edge, Che a8Susq>clon that pressure is constant across the bourxiary layer Is no 
longer valid. Thus, the elimination of these singularities should allow the 
viscous flow on an airfoil to be owre reliably predicted using conventional 
boundary layer theory* More significantly, because of the reduction in the 
strong viscous effects at the trailing edge. It Is reasonable to expect that 
airfoil aerodynamics might he enhanced as the level of performence ach.eveu In 
the real flow field would more closely approach the high levels predicted 
using potentlel flow methods* Finally, because Che msnner In which the flow 
approaches the trailing edge should be globally influenced by the removal of 
the trailing edge singularities, It Is anticipated that the possibility ot 
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large upper and lower surface velocity differentials in the vicinity of the 
trailing edge would be eliminated. Thus, separation problems introduced by 
the steep adverse pressure gradient that typically exists on the upper surface 
near the trailing edge when such velocity differentials exist would be 
mitigated. 

In order to examine the flow in the vicinity of the trailing edge in 
detail, consider the transformation of the unit circle centered at the origin 
of the emplane into an airfoil in the z-plane as shown in Figure 21. The 
complex potential function for the unit circle having circulation T and an 
angle of attack a is 

F(5) - U - U ^ In ; (3-1) 


and the complex velocity is 

- U - U + 'll (3-2) 

Imposing the Kutta condition, the circulation necessary to fix the stagnation 
point at ^ i in the circle plane is found to be 


r ■ 4 tt U sin a 


(3-3) 


Thus y 



glo^-2 ^ 215 “^ sin o 


U 


(3-A) 
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which may be written 


w(^) 

U 


-ia 


.‘ V ‘)(1 - 


r‘) 


(3-5) 


The velocity in the airfoil plane is 

w(z) = u-iv * ^ “ v(0 / z'(c) 


(3-6) 


To preclude the presence Of a stagnation point at the trailing edge of the 
airfoil, Zq., it fs necessary that the trailing edge of the airfoil be 
cusped. This requires that 

2^(C) » 0 (3-7) 

Because the velocity at the trailing edge, from equation (3-6), is of an 
indeterminant form, L'Hopital's rule can be invoked to yield 


w(z^) 



(3-8) 


At this point, the complex pressure gradient is defined as 


R 





(3-9) 


Using the Bernoulli equation, the pressure at a point on an airfoil, p, may be 
related to the free-stream conditions such that 


p • P [U^ - (u^ + v^)] 


(3-10) 
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(rfhich is dif ferenclaced to give 

= -p (uu^ + vv^) (3-11) 

Py * “P [uUy + vVy) (3-12) 

By substituting equations (3-li) and (3-12) into equation (3-9), it is found 
that 


R - -P [u(u^ - iUy) + v(v^ - iVy)] (3-13) 

Continuity and Irrotatlonality require that « Uy and “x * “ ''y' Thus, 
equation (3-13) nay be written 

R « -p [u(u - iv ) + vfu ~ iv 1] 

I V X x' ^ y y'-' 

* -p [uw' + v(iw')] 

■ -p [u + iv] w' (3-14) 

or, more simply 

R ■ -p w(z) w'(z) (3-15) 

where w(z) is the complex conjugate of w(z). At Che trailing edge, the 
complex pressure gradient is given by 




R<P ■ — p w( Zip ] w ( z j ) 


(3-16) 
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In order to evaluate equation (3-16), an expression for w'(z^) is 
required. Note that the result of differentiating equ ition (3-8), 








(3-17) 


C ♦ C, 


is of an indeterminant form such that L'Hopital's rule is applicable. Because 
w(^_) and Zj' are both zero, this yields 


w'fz^) « 11m - y(z)z . ”* . 

T C > 3rz-)^z" 


(3-18) 


which is still of an indeterminant form. Using L'Hopital's rule again gives 


w’(zt) “ 


rw'*(e)z‘‘ - w*(c)z ' ' ' 


? - Cj. 6(z'](z")^ + 3(z')^z'" 


(3-19> 


Now« because the denominator of this expression is zero, v'(z^) will be 
unbounded unless the numerator is also zero. Thus, for w'(z<^) to be finite, 
it is necessary that 




If this condition is satisfied, then equation (3-19) will be of an 
indeCerminant form and the use of L'Hopital's rule again yields 


(3-20) 


♦ 

1) ■ 

I i ' 

I 

i ; 

iii 




w'"(Cr)zf - w'(c^)z., 


(3-21) 


Now, subscicucing the above relation, along with the conjugate of equation 


% 
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(3-8), into equation (3-16), the value of the complex trailing edge pressure 
gradient, provided that equation (3-20) is satisfied, is given by 


“ “P 




3 


y [zf’ (3-22) 


To further evaluate equations (3-20) and (3-22), equation (3-5) may be 
differentiated successively and evaluated at the trailing edge giving 


- = 2cosa (3-23) 

— 0 * - 6cosa - 2isina (3-24) 

■ 24cosa + 12lslna (3-25) 


By writinj equation (3-20) as 

Zt"’ _ 
Zt" " 


(3-26) 


and making use of equations (3-23) and (3-24), the condition required for the 
complex trailing edge pressure gradient to be finite on an airfoil may be 
expressed as a condition on the transformation function requiring that 



3 - itana 


(3-27) 


While the ramifications of this requirement will be considered later, note 
that the iisaginary part can only be satisfied for a particular airfoil at only 


i 
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a single angle of attack. At any other angle of attack, the complex trailing 
edge pressure gradient is unbounded. 

Using equations (3-23) and (3-25) in equation (3-22), the value of the 
complex trailing edge pressure gradient is given by 

R a - i ?L [6z^(2cosa + isina) - z ^'^cosaj (3-28) 

" 3 (^ 7 ) ^ 

provided that the condition of equation (3-27) is satisfied. Non- 
dimensional izing equation (3-28), where £ is the chord length of the airfoil, 
yields 


■ SxT - ‘ Vt 

» ^ [6z_“(2cosa + isina) - z_^'^coso] (3-29) 

3(z^")(z^")^ 

By considering a circle of radius r and center location u in the C“plane, 
rather chan a unit circle centered at the origin, the results of this analysis 
can be made applicable to the von Mises airfoils previously considered. In 
that case, it is interesting to note that all of the symmetrical airfoils 
generated using the minimum circle radius which generates a physically 
realizable airfoil, satisfy the condition for a finite complex trailing 

edge pressure gradient. In order to explain this, consider that, in the case 
of overlapping airfoils, the point of cross-over moves aft as the generating 
circle radius, r, is increased. When r • ^^niin» cross-over point is at the 
trailing edge. Thus, the trailing edge of the airfoil is locally like the 
flat plate at zero angle of attack and, evidently, the flow behavior is 
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similar. As an example of such an airfoil, the von Mises profile depicted in 
Figure 22 has trailing edge pressure gradient values of 

Cp^.p » -1.37, while, as is the case for all non-lifting symmetrical airfoils, 
CpYT 0. For comparison, the airfoil of Figure 16 has Cp^-j* * -.12 and CpY*j ■ 
0. It should be noted that because having r * case of two 

generating zeros yields the Joukowsky flat plate. Figure 13, that this is the 
only Joukowsky airfoil which can satisfy the condition for bounded trailing 
edge pressure gradients. This fact explains the presence of the "hooked" 
shape observed at the trailing edge of Joukowsky airfoil velocity 
distributions. 

The Eppler Airfoil Design Method 

In order for airfoils to be developed which make use of the conditon for 
finite trailing edge pressure gradients, it is necessary to incorporate the 
necessary condition into some airfoil design methodology. To this end, the 
inverse scheme developed by Eppler was selected as being the most suitable for 
this purpose. The theoretical details of this design procedure are presented 
in Eppler [5], and also summarized in References (8], [17], and [19], while a 
description of the code and its usage are documented by Eppler and Somers 
[23]. As modifications required to allow the method to design airfoils having 
finite trailing edge pressure gradients are significant, before considering 
these changes, a summary of the method will also be included here. 

As depicted in Figure 21, the Eppler airfoil design procedure Is based on 
the conformal mapping of the flow field exterior to a unit circle in the C- 
plane into the flow field exterior to an airfoil in the z-plane. The complex 
velocity for the flow about the unit circle is given by equation (3-5) where. 
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to be consistent with the derivation of Eppler, the free-stream velocity is 
taken to be unity. The mapping function, z(c), naust, as detailed in the 
discussion of the von Mlses transformation, preserve the flow conditions at 
infinity. Thus, it is required that 

z(®) - » and ^ (3-30;, 

Most generally, these requirements are met by the power series 

00 

zU) - BjC + I * 0, real) (3-31) 

v-0 

where the 6 *s are limited to values for which the series is convergent when 

V 

|c| > 1. As consistent with equation (3-6), the velocity In the airfoil plane 
Is given by 


w(z) 


-19 dF dF dc 

dz dc dz 


(3-32) 


This relationship, for reasons discussed by Lighthill [9], is more 
conveniently represented by 


In w(z) ■ In 


ie 




(3-33) 


Because In -rr is known from equation (3-5), chen if v(t) is specified, the 
dC 

d dF 

real part of In ^ follows directly* Furthermore, because In ^ is analytic 

In the exterior of the unit circle then, using the real part which is known <>n 

the boundary, C “ e^^, the imaginary part can be determined. Thus, the 
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function — , which maps the circle into the airfoil, can be solved for from 
dC 

equation (3-33) and z obtained by integration. The advances which are 
represented by the Eppler method include the development of an algorithm for 
the numerical computations required by the procedure just outlined, and the 
Introduction of a form for v(4>) which allows practical airfoils to be obtained 
Ln a straightforward manner. 

The numerical algorithm of the method is based on the Introduction of 
dz 

In — in the form given by 

i (. ♦ibjr" (3-34) 

m-0 

This form is advantageous in that it isolates the singularity which occurs on 
the boundary, ^ * e^^. Operating on equation (3-34) yields 


dz 

dC 



l (* 


m*0 


lb.)4 


-m 


(3-35) 


This result must be consistent with equation (3-31). Thus, differentiating 
equation (3-31) and comparing terms having like*powers of t to those of the 
preceding expression, it is found that compatibility requires 


*>0 

- 0 

(3-36) 

*l 

- 1 

(3-37) 


- 0 

(3-38) 


Equation (3-36) is necessary in order that there is no rotation of the 
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free-stream velocity at Infinity* Equations (3-37) and (3-38) make 
the C ^ term vanish and, as a consequence, ensure that the profile will close 
at the trailing edge* Additionally, the requirement that velocity of unity be 
preserved at infinity, as given by equation (3-30), necessitates that 

- 0 (3-39) 

Evaluating equation (3-33) on the circle boundary, t and rearranging 

the result gives 

[in ^ - In f 1 - i)] - In v(4) + 16 

04 4 C-e ^ 

» in - ln(l - 

- P(^) + i Q(^) (3-40) 

P(i^) and Q((^) are defined as Che real and Imaginary parts, respectively, of 
this expression. The use of equations (3-5) and (3-34) allows P(<>) to be 
expressed as 

P(Q) ■ I (a cos + b sin m^) • - ln[ 1 (3-41) 

B-0 " B ^ 2|cos(|-o)| 

so Chat Che coefficients, Che a^^'s and b„'s, are determinable In the same 
manner as chose of a real Fourier series. Consequently, using Che formulas 
for Che evaluation of Fourier coefficients, and noting Chat Che term 
containing bQ disappears In the expansion of equation (3-41), the requirements 
of equation (3-37)-(3-39) generate three integral constraints. 
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7T 


ir 


^0 


27f 

( 

J 

0 


27T 


/ 

0 

2tt 

/ 

0 


P(^) d(^ » 0 
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?(<>) cos (t “ n 


P('p) sin di^ * 0 


(3-42) 

(3-43) 

(3-44) 


which nusc be satisfied by Che specified velocity distribution, v(t). In the 
actual loeChod, it is not necessary to determine the rem/iining series 
coefficients but, rather, necessary only to calculate the conjugate harmonic 
funcion, Q(i^). Thus, using Poisson's formula, Q(<^) is given by 


Q(*) 


I 

m*0 


fb cosm^ ~ a sinm^) 
n ni 


2it 


/ P(>|/)ctn d<» 

0 ^ 


(3-45) 


Substituting equation (3-41) into equation (3-40) and simplifying yields 




i(Q - |) 


(3-46) 


and, for dC • ie^*^ d^. 


i£.n riii 


-4sin I |cos(|- - o)| e 


l(Q ♦ I) 


(3-47) 


Splitting this expression into its real and imaginary parts gives 
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^ - Asin I |cos(| - a) I cos[| + Q(«)] (3-48) 

^ - Asin I |cos(|- - a) I sin[-| + Q(^)] (3-49) 

Thus, once Q('>) has been determined, the airfoil coordinates can be obtained 
by simple quadrature. 

Up to now, the derivation has proceeded In such a way chat a fixed angle 
of attack, a, is selected, as well as a velocity distrlburlcn, v(t), which 
must satisfy the three Integral constraints, equations (.3-42)- (3-44). Then, 
using equations (3-48) and (3-49), a profile is determined whose velocity 
distribution, v(<^, o), agrees with the prescribed v(<1>) at the selected angle 
of attack. The method, however, is considerably more flexible. From equation 
(3-35), it is clear that Che transformation, and therefore the airfoil shape, 
is fixed once Che a„|'s and b„'s have been determined. As a change in a has no 
effect on Che coefficients, P((^) is independent of a. Consequently, equation 
(3-41) can be written 


P(9) 


- In ^ 


2 1 cos 


v*(^) 

(f ■ “*) 


(3-5U) 


where v*(^) is the velocity specified at a point on the profile corresponding 
to end a*(4) is the angle of attack at which that v*(i>) is to be 
realised. Thus, Che velocity distribution v(^, a), for any arbitrary a, is 
obtained from equation (3-50) an 


I cos (y - «)| I cos f| - a*) I 


(3-51) 
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where is invariant with respect to a. While v(<^, a) must be continuous 

over the a^’^foil, it is possible for a discontinuity in a* to he compensated 
by a discontinuity in v*(t)* Thus, rather than specifying only a single a* at 
which the velocity distribution v^(<>) occurs, unlike other inverse methods, it 
is possible to select different values of a*, as indicated by the notation 

for different segments of the airfoil. Consequently, different parts 
of the airfoil can be designed for different angles of attack. In this 
manner, an airfoil can be designed from the onset to have the desired 
aerodynamic characteristics over a chosen range of flight conditions, rather 
than being point-designed for a single angle of attack and mod<>fying the 
result until acceptable off-design performance is achieved. 

As mentioned, the numerical method is dependent on the convergence of the 
series contained in equation (3-34). This convergence is assured if the 
velocity distribution, is specified such that it makes P(4)f equation 

(3-41), a piecewise continuous function containing, at most, points having 
finite jumps in the first derivative* In this case, as detailed in Reference 
(28], the a^i's and b^’s are 0(l/m^)* This not only allows successful 
mimtrical treatment, but also guarantees that the resulting profiles are 
smooth* 

In order to implement the specification of the velocity distribution, the 
unit circle in the t^plane is divided into Ip segments over the Interval (0, 
2w), Referring to Figure 23, 

♦o ■ 0 < ♦! < ♦a ^ ^ ♦! <3-52) 

P 

and, in addition, denotes the arc limit containing the leading edge 
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stagnation point achieve the flexibility Introduced by being able to 

specify different values of a* for different segments of the airfoil, 
o*(^) will be considered as constant over each of the defined Intervals such 
that 


a*((^) * “ constant '^i-1 ^ ^ ^ (3-53) 

The specification of the velocity v*(<^) takes the form 

v^(^) - v^w(^) (3-54) 

where v^ Is taken ae constant within each circle segment, and 

w('>) is a function which Includes a term allowing for main pressure .ecovery 

on the airfoil, as well as a contribution to Insure that the profile closes. 

On the upper surface 0 < , the form of w{<^) Is 

1 , A e 

cos^ - coa^ -u cos^ - cos<> 2 H 

'<♦' • t' * rv c..» “ 1 1 - •“< I - co.« 

w ’s 

while on Che lower surfece, < 2e, chi peraneCers u, , and 

> are replaced »y u, t > end ^ , reapecClvcly* The expressions within 

s n w 8 

Che braceSi {f(d)}i ere creeced ea special functions in chati if f(>) < U then 
(f(^)} • 0 and if f(^) > 0, {f(d)} ■ f(^). For siaplicicy. equation (3-55) 
for the upper surface oiay be written 

K 

w(i) - W^C*)'" W^(4) " (3-56) 


Mi 


w# a 
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rind, for the lower surface. 


w(^) = W^(^) W (^) ” (3-57) 

w s 

In both cases, the W (<j>) term produces the main pressure recovery, the length 

w 

of which is governed by the specification of 4> . The total amount of recovery 

w 

and its particular shape are established as part of the velocity 
specification, along with the and ot^*s, by the parameters K and p. In 

the term generating closure of the airfoil, the quantity is fixed by the 

specification of while the exponents, Kj| and are left free to be 
determined by the solution procedure, A typical velocity distribution, 
-^(<)>(x)), is sketched in Figure 24* 

Substitution of equation (3-54) into equation (3-50) yields 

P(<^) * Injcosfl'- ] I - In v^-ln w(4>) + ln2 (3-58) 

(4>J_1 < t ; 1=1, 2, 3 Ip) 

Thus, ac the trailing edge, continuity of P(^) requires that 


InjcosaJ - In v^- In w(0) • Ln|cosOj |“lnvj - In w(2n) 

P P 


(3-59) 


whereas at all other segment boundaries. 


\ ^1 
ln|cos(2 Oj)|-ln - In w(ifr) - ln|cos('^ ®1+1 
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- In - In w((Ji) (3-60) 
; i-1,2.3 Ip_^) 

A problem arises with the velocity distribution 
stagnation point, given by + 2a^, in that P(<J>) 
point and no longer satisfies the conditions imposed 
be avoided, however, by requiring - 2a^ le ^ ^ 
or, more simply, by requiring 


at the leading edge 

is undefined at that 

on it. This problem can 

and , - 2at. , 

i,le i,le+l ’ 




i.le 


(3-61) 


Ac this point, Che method requires chat Che values of Che pressure 

recovery and closure parameters, u, u, i , > , and ()> , be given, along 

w w s s 

with all of the pairs. The Ip constants, the v^'s, and the closure 

contribution exponents, K|| and K^, are solved for such that the Ip matching 
conditions, equations (3-59) and (3-60), as well as the three integral 
constraints, equations (3-42)-(3-44 ), are satisfied. Consequently, as given, 
the problem is over-specified and it is necessary to relax one of the given 
parameters so that all of the required conditions can be met. Because of its 
strong influence as the matching point of the upper and lower surface velocity 
distributions, the method uses the leading edge arc limit, as the 

necessary free variable. 

By substituting equation (3-58), along with equation (3-56) or (3-57), 
into equations (3-42)-(3-44), it is found that the integral constraints can be 
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evaluated in closed form. With the substitutions indicated, equation (3-43) 
becomes 


2it 2it 

/ P(i|»)cos(|) d(j) = / [ln|cos(-j - a.]| - In v.(i))) 

0 0 ^ ^ ^ 


+ U In In + In 2] cos<J> di|) = tr (3-62) 


Now, defining 


2tt 

- / In W (4>) COS4. d(J> * W , (3-63) 

0 s Cl 

2it _ 

- / In W (i) cosiji diji =• W , (3-64) 

0 ® p 


and, introducing the notation 


ln(i,j) » In I cos (-^ - Oj)| (3-65) 


The evaluation of equation (3-62) gives 

2it _ 

f P(<p) cos<^ dvt- - + Kjj W J + 1 {sin 2a^ 

0 p i-l 

lln(l,l) - ln(i-l,i)J 


+ j (sin - sin ) + sin (>^[ln(i,l) - In v^j 
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- sin 4 )._ Jln(i-l,i) - In v.]} 


_ 2n _ 

+ u / In W (<p) co8i> + V j In W (<>) cosi> d4 * n 

0 ” ? ” 

w 


(3-66) 


Using the matching condition, equation (3-60), it is found that a number of 
terms in this expression drop out such that, after simplification, it may be 
written 


S. “cl * "cl * -"c ■ 0 
P 


(3-67) 


where is defined as 


•J- * I {sin 2 a [ln(i,D - ln(l-l,i)] 

c £.1 1 


1 ” 

T ” '^1 I )co® 2o. } - X + w / In W (<^)cosit> 
^ 1 11 w 


d4> 


2x 


+ U / In W„(ifr) cos^ di(> 


(3-68) 


Now, after formally setting Oq and a^ to zero, the terms containing iji^ are 

P 

collected and the leading edge arc limit, is isolated such that 

equation (3-68) is rewritten as 


J ■ a + b ln(le,le) + c In (le, le + 1) + d , (3-69) 

ccc c citie 


where the coefficients are given by 
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I , 1 e- I 

a = {sin 2a. ln(i,i) - sin 2o. . ln(i,i+l) 

i-0 ^ ^ 


+ y t. (cos 2a - cos 2o , ]| + \ 


]| + I {sin 2a. ln(i,i) 

i=l ,le+l 


“ sin J (cos 2a^ - cos 


. 2tt _ 

- 7T + Vi f In W ((Ji) cos^ d<() + y / In W (((>)cos4> 
0 ^ 

w 


b = sin 2a. 
c le 


* “ sin 2a, 


T (cos 2a, - cos 2o, 


Using Che definitions given by 


/ * In W (i>) sini^ di> ■ W 
0 ® * 


/_ In Wg(li) sin(|) d4i - W^j 
♦g P 


the Integral constraint of equation (3-44), in a manner similar to that 
applied Co equation (3-43), becomes 


*n “si * 5. “si * J. ■ 0 
p 
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where Jg is ultimately expressed as 


•^s ” + c ln(le,le + 1) + d 


s ^le 


(3-77) 


and the coefficients are 


i,le-l 


a = 
s 


I { -(l + cos 2a ) ln(i,i) + [l + cos 2a,^,l 
i»U ’■ i+l'^ 

ln(i,i+l) + j (sin 2a^ - sin 

I 

P 

L { “(^ cos 2a ]ln(i,i) + (l + cos 2a , ) 
i,le+l ^ 

ln(i,i+l) + j (sin 2a^ - sin 

” _ 2it _ 

+ U / In W^((ji) sin0 d(fi + u f_ In W ((J>) sini> d(|) 

0 (fc 

w 

b * - ( 1 + cos 2o, ] 

s'- le-* 


(3-78) 


(3-79) 


• (l + cos 


(3-80) 


-j (sin 2o^^ - sin 2o^ ) 


(3-81) 


At this point, consider the trailing edge matching condition given by 

equation (3-59). By repeated applications of equation (3-60), the velocity 

terms v^ and Vj can be eliminated and the resulting expression written as 
P 
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- Kj^ In W^((J) + Kjj in W^(2 ti) + J^. = 0 


(3-82) 


in which 


ln(le,le) + ln(le,le+l) 


(3-83) 


and 


i,le-l 

a ^ {- ln(i,l) + ln(i,i+l)} 

^ i«0 


I 

P 

* I (- ln(i,l) + ln(l,i+l)| - u In W (0) 

i-l,le+l 


M In W (2 it) 
w 


(3-84) 


\ 


(3-85) 




(3-86) 


Thus, equations (3-67), (3-76), and (3-82), can be used to solve for the 
three unknowns, K|j, K^, and 6^ In matrix notation, this system of 

equations may be represented 
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W 



(3-87) 



Dj - W^j In W^(2ir) + In W^(0) (3-93) 


D, - - W^, In W^(2ir) - W . In W„(0) (3-94) 
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Equation (3-88) can now be solved numerically for using the 

technique of Regula falsi. Once has been obtained, Jg, and can 

be calculated and used in a straightforward manner to solve for and K^. In 
order to determine the v^*s, the first Integral constraint, equation (3-42), 
is evoked along with the matching conditions given by equation (3-60). Thus, 

2tt ^l V, 

I P(4>) d4> - I 1 / [ln|cos {j - a^^)! - In 
0 i»l I 

+ u In W (♦) - K„ In W (4>)1 d(ji + 2n(ln 2 - In v, ) (3-96) 

W n S 1 

The velocity ratios, — , are found from sequential applications of equation 

''l 

(3-60) and the results used to determine v^ from equation (3-96). At this 
point, P(i^) and Q(6) can be determined and the coordinates of the profile 
follow directly. 

The closure requirements of the airfoil determine the values of Kj| and 
which, in turn, influence the trailing edge closure angle. In order that 
some control over the trailing edge shape is possible, the digital computer 
coding of the method allows a value of K^, where * K|| -»■ K^, to be 
specified. Then, any one of a number of iteration schemes can be selected 
which vary particular combinations of the specified Input parameters until the 
desired value of is achieved. This procedure, along with a detailed 
discussion of the numerical Implementation of the method. Is contained in 
Reference (25]. 

Incorporation of the Conditions for Finite Trailing Bdge 
Pressure Gradients into the Eppler Design Method 


As developed previously, the condition on the transformation function 
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which results In a finite complex pressure gradient at the trailing edge of an 
airfoil is given by equation (3-27). In order to Incorporate this result into 
the design method of Eppler, the representation of the derivative of the 
transformation function, provided by equation (3-35), is differentiated to 
yield 


m-0 


-m 


z ■ e 






(3-97) 


which, when evaluated at the trailing edge, C ■ 1, becomes 


I {a ^ ib 1 
m*0 

■ e 


(3-98) 


. i 


Differentiating equation (3-97) yields 


I 

»■ 1 


V (a ib )c"“ 

I ^ m m-' 


m“0 


z “ e 


' ' 




m»0 


li ; 


+ (I - t) 1 ^ 

ofO 


-(m+2) 


’ i ( 






(3-99) 


For this expression to be evaluated at the trailing edge. It is necessary that 
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the last two terms drop out as occurs if the summation present in each term is 
finite. Noting that the and can be considered as the coefficients 

of the Fourier series representing P(<^), equation (3-41), and, considering the 
theorems concerning the differentiation and convergence of such a series. 
Reference (28j for example, it is found that the summations will be finite if 
P’(<t) is continuous at the trailing edge* Thus, ii is required that 


P’(0) - P*(2ir) 


(3-iOO) 


If this is the case, then equation (3-99) becomes 


m»0 


2 - 2 I ra(a + ib )} 
^ A ro 

m-0 


(3-101) 


Substitution of equations (3-98) and (3-101) into the condition for finite 
trailing edge pressure gradients, equation (3-27), yields 


' ' ' m m 

^ — • -2 - 2 \ ma - 2i \ mb ■ -3 - itana 

T m*0 ™ m*0 * 


(3-102) 


To further resolve this expression, the series representations of P(^) and 
Q(^), given by equations (3-41) and (3-45), are differentiated to obtain 


1 


f(*) - I 

m"0 


- aa_ sin oi4 mb. cos 


(3-103) 


Q'(4) - I 

a»0 


I cos w4 ■“ mb. sin a4 

■ n 


(3-104) 


't) 


4 • 


* a 
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which, when evaluated at the trailing edge, become 


P'(0) - mb (3-105) 

m*0 

m 

Q'(0) - y ma (3-106) 

m*0 

Using rhese results, equation (3-102) yields the conditions that 


P'(0) - tan a 


(3-107) 


Q’(0) - - i 


(3-108) 


which nust be satisfied in order for a airfoil to have finite trailing edge 
pressure gradients. 

In order to consider equation (3-l07> further, the expression for P(<>) of 
equation (3-58) is differentiated and evaluated at the trailing edge yielding 

f(0> (3-10») 

In asking Che transition froa v(<^) to V(x), it is found that the velocity 
distribution on the airfoil, V(x), has an infinite slope at the trailing edge 
unless w'(0) ■ 0. Because of this, it already a requireaent in the aethod 
that w'(0) ■ 0 and, consequently, equation (3-109) becomes 

• -y can Oj 


P*(0) 


(3-110) 



Equating this result with equation (3-107), it is found that the condition 
imposed on P'(0^ is satisfied when 
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a - Oj (3-111) 

Thus, the complex trailing edge pressure giadient on an airfoil can only be 
finite when the airfoil is operating at an angle of attack corresponding to 
that specified in the design process for the first arc segment. 

To satisfy the second condition required for bounded trailing edge 
pressure gradients, equation (3-108), the integral expression for Q(t), 
equation (3-45), is differentiated and calculated at the trailing edge giving 

, 2x 

Q’(0)-^/ P(tJ-) csc^ (f] d* (3-U2) 

^ 0 ^ 

Substituting this result into equation (>-108) introduces an additional 
Integral constraint which is given by 

2 « , , 

/ P«r) CSC f d» - - 2t (3-113) 

0 

In addition to the conditiona of equations (3-111) and (3-113), 
examination of equation (3-22) reveals that in order for to be finite, the 

4 u 

fourth derivative of the tranaformation evaluated at the trailing edge, , 
must alao be finite. Evoking the same arguments used for insuring the 
flnlteness of it found that mill be bounded provided that P“(|) 

Is continuous at the trailing edge. To insure this, consider the result of 
differentiating equation O-M) twice, giving 
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_ I ci _ w”(<» + 



0-114) 


[«i-i < ♦ < ; 


1 - 1 , 2 , 3 ,. 


1 

P 


1 


Because w'(0) - w'(2n) - 0, It Is evident that the continuity of P”(<>) at the 
trailing edge requires 


w-(0) . w-(2ir) 
w(0) w(2w) 


(3-115) 


After perforaing the necessary manipulations, making use of equation (3-53), 
it is found that this condition can be written more specifically as 

uK . , ^ 

1 + K (l-K)cos ♦ M - cos ♦ ^ 

W 8 

- V 

• JiLJs _ ♦ 1.125 ( 2-^) (3-116) 

1 + K ♦ (l-K)eos 6^ 1 - cos 6^ 

It should be noted that satisfying this condition also guarantees that the 
last two terms in equation (3-99) drop out. 

As presented thus far, the development has shown that an airfoil which 
satisfies equations (3-111), (3-113), and (3-116) will have finite trailing 
edge pressure gradicncs; however, in order to incorporate these conditions 
Into the design process it remains to evaluate the integral constraint, given 
by equation (3-113). in terms of the appropriate design parameters* To begin, 
although the integrand of equation (3-113) is singular, the Integral can be 
shown to exist in the Cauchy Principal Value sense by d If tsrent fating the 
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expression for y(f)» equation (3-45), to obtain 


2n 


Q'('i') * / P('l') csc^ di|» 


^nd, integrating oy parts, yields 


Q’(t) - - P(2ir) ctn (it - |) + ^ P(0) ctn (- i) 


2^ 2n 


2it 


+ -^ / P'('l») ctn dt|» 


0 


Because the required continuity of P(^) gives P(0) ■ P(2v), the 
of this expression cancel leaving 


2ir 


Q'(4») “ / P'(’l') ctn (-^^^) di(> 


and Che fou :h ntegral constraint Chen becomes 


2x 


Q’(0)-^/^ P'(>l<) ctn I d* • - j 


After renaming Che variable of integration, this result is used 
as 


2x 


I, - 2 / P'(^) ctn 4 d* - - 2W 

^0 ^ 


(3-117) 

(3-118) 
leading terms 

(3-119) 

(3-120) 
Co define I 4 

(3-121) 


Now, differentiating equation (3-S8) '[Ives 
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P’(«) 



w'(i|i) 

w(({i) 


!*i-i ‘ 


0 < ; i-l,2,3,...,Ip) 


(3-122) 


Substituting this expression Into equation (3-121) ultimately yields a result 
which Is denoted as 


where 


1,-1 + I. - - 2x 

4 a b 


(3-123) 


~ 1 I tan (-j - a ) ctn 4 d<t> 

1-1 * 1-1 


(3-124) 


and 


; 1 - 1 , 2, 3 1 1 






(3-123) 


The analytical evaluation of is now undertaken , beginning with 
consideration of the indefinite integral 


laj ■ / tan - o^) ctn d<^ 


(3-126) 


for which, after some aanlpulatlon. It la found that 


Nil « 
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I , * - 2 tana. In 

al 1 


sin j 

■li - »i) 


(3-127) 


Using this expression in equation (3-124) gives 


^ “ 1 ^ 2tana. In | 

^ 1 _ 1 1 


sin j- 


i=l cos(y - a^) ‘^^i-l 


(3-128) 


which becomes 


P y 

^a * ^ {■ •I'i 't'j.j + 2tana^ [injsin 


- ln|sin -^i| ] 


T • ™ 1 

- 2tana^ [ln|cos - a^)l - ln(cos - a^)|j} (3-129) 


Performing the Indicated summation with * e -► 0 and * 2ir yields, 

P 


9 

Ig * 2tano^ [ln|sin - Injsin yj ] 


1>i 

- 2tana^ [injcos “l ^ “ ln|cos(-a^ 


2w + 2tan Oj [ln|8in(w " j] ” Injsin 
P 


■0 


- 2 cana^ [lnjcos(u - )j - Injcos ( )j 

P P p 

P ^ J ^ 

+ I ( 2 tana [injsin - 5—1 - Inj sin — ^~ -j J 
i -2 ^ ^ 


J 


. .4 
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Extracting the singular terms in the above expression and combining them as 
Sxa gives 


Sj » - 2tana^ In|sin y| + 2tana^ Injsin fx * y] | (3-131) 

Because |sin(x "" y)| * I sin y| , the singular terms can be eliminated by 
requiring that 


«X * (3-132) 

P 


The expression for that remains is then simplified by incorporating 
additional terms under the summation such that 


p — 1 4^ 

I « - 2it ♦ 'I {2cana. [ln| sin - Injcos (- 5 ^ - 0^)1] 
a 1 ^ z 1 


T 9 

+ 2 tano^^j [ln|cos - Oi+i)l " ln|sin-^ 


i +1 


(3-133) 


To begin the evaluation of Ijj, equation (3-125), recall that w'((>)/w((>) 
varies In a piecewise fashion with Using the form of equation (3-36) to 
describe w( 6 )> and comparing that with equation (3-55), it Is found that, on 
the upper surface. 


W„(6) - 1 + K (-j-^ 


C 08 ^ - C 08 ^ 


C 08 ^ 




(3-134) 



r 
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Wg(4>) = 1 - m (-p^ 


cos4> - cos<|) 2 


cos^ 


(3-135) 


while, on the lower surface, 


_ _ cos<> - cos ^ 

w^(^) = 1 + K [ 

1 + cos ^ 


(3-136) 


_ COS^ - cos 

w^(4i) ■ 1 - m ( — -] 

1 - cos ♦ 


(3-137) 


With regard to Che earlier discussion concerning the piecewise creatmenc given 
to the function described by equation (3-55), Ijj Is evaluated over the 
piecewise segments defined as follows: 


i: 


W*(<») 

v('» 


w'(») 

w(^) 


w •(♦) W •(<)) 

wT^ wYTT" 


- u 


Ww’( 4 >) 

«„(♦) 


0 < 4) < 


L < ♦ < * 
s w 


(3-138) 


(3-139) 


I 

0 

D 


w( t>) 




w(!>) 


V*) 


w’(6) 

2L^ • - y + V -2 

"W(^) 


«.<♦) 


() < 4) < A 

w ^ ^w 


^ ^ <b ^ A 

W ^ ^8 


^ t < 2w 
$ 


(3-140) 


(3-141) 


(3-142) 




^ • 
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Thus, for Che mosc general case, Che Indeflnlce Incegral Is given by 


where 


and 


lu " Iwi + luo 
b bl b2 


( 3 - 143 ) 


‘bl ■ ^ in?r "" f 


( 3 - 144 ) 


ct) 


w '(♦) . 

^b 2 * ~ ^ ^ W ( 41 ) 2 

s 


( 3 - 145 ) 


In Che following evaluacion of only che noCaCion for Che upper 
surface will be used; however, Che resulcs are valid for Che lower surface as 
well, if the upper surface quantities are replaced by their appropriate lower 
surface counterparts. In order to evaluate let 


* « 


® 1 + cos^ 


( 3 - 146 ) 


so that 


f - 1 - g cos ♦ 


( 3 - 147 ) 


» T 

S 


*„(♦) * f + g COS(> 


Thus, the integral becomes 


( 3 - 148 ) 


» 3 
. i 


1 

I • 


I m - 2u I f— a-2isi_1 ctn 4 d« 

bl ^ + g cos^^ I ’ 


( 3 - 149 ) 


1 i 




'I ' 

» \ 
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After some manipulation, it is found that 


- 2u<> + UF [<P, f, g) 


(3-150) 


where 


2(s-n ml 


[g-f Jtan -J + ig - f J 


(T-. *_(777i^ 


(3-151) 


-1 I , 

(^2.T)T72“» . ,2, -.Til 




*2.2 

f > g 


(3-152) 


In order to evaluate Ijj 2 » l-®t 


N -_iS! (3-153) 

1 1 - cos ♦g 

a - 1 - /« (3-154) 

b * 1 + /m (3-155) 

so that 


Wg(4) * [a + (1 - cos4)][b - (1 - cos^)] (3-156) 

and 

W^'(4) • Nj 8ln6(b - Nj(l - cos*)] - sin*(a + Nj(l - cos*)] 

(3-157) 
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The expression for Ijj 2 nay then be written as 


^b2 * ^b2a * ^b2b 


(3-158) 


where 


^b2a “ " ^ [a + Nj(l - Ls(|i)] 2 


± 


di|i 


(3-159) 


b2b 


^ [b - cos*)] 2 


1 


d()i 


(3-160) 


After performing the integration indicated in equation (3-159), ljj 23 becomes 


a + 2N, 1/2 , a + 2N, 1/2 . 

^2. ' V - “Si t1 


(3-161) 


Evaluation of the expression for Ij,2b» equation (3-160), after some 
manipulation, gives 


‘b2b ■ * “Si V*- "i- '>> 


(3-162) 


where 


(2Nj - b) (b - 2Njtan | 

h “ 7:1 ~rr72 


(b - 2Njb) 


(b - 2Njb) 


rrr 


(3-163) 


(2Nj - b) 


In 


(2Nj - b)tan + (2Njb - b ) 


2U/2 


2(2bNj - '(2Nj - bjtan - f2Njb - b‘) 


irrTil ^ ^ 


(3-164) 
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Combining equations (3-150), (3-158), (3-161), and (3-162) Into the 
expression for Ijj, equation (3-143), yields 

« - 2u^ + U F(<>,f,g) + 4 Kjj(> + 4K^ Kp((fr,N^,b) 

a + 2N, 1/2 , a + 2N, 1/2 

■ ^ ^ i ^ (3-165) 

Using this result, can be evaluated at the appropriate integration limits 

for the given Intervals so that 

a + 2N, 1/2 , a + 2N, 1 fl 

'b ■ % % (-T— 5^1 

+ Kp(<>^,Nj,b) - 2p0^ + uF[^^.f,g) + 2 u 

- iT F f, i) - Alt u + 8 it - 4 - 4 Nj. b) 

_ a + 2N, 1/2 , a + 2N, 1/2 T 

+ 4Kjj (— ^ i-) tan‘M(-^ tan j^] (3-166) 

This result, combined with that obtained for equation (3-133), provides an 
analytical representation of the Integral constraint, I4, which must be 
satisfied In order to Insure that the trailing edge pressure gradients of an 
airfoil are bounded. 

At this point, there are four Integral constraints and Ip matching 
conditions which must be solved for K)|, K^, and the Ip unknown 

constants, the Vj^'s. Thus, In order to satisfy a fourth constraint, an 
additional parameter must be relaxed. Because for most applications of 
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Interest, the specification of the upper surface velocity distribution is of 
greater consequence than the specification of the lower, both u, and K, which 
control the shape of the lower surface recovery, are candidates; however, as 
freeing u results in a more direct method of solution, it is chosen as the 
additional parameter to relax* 

At this point, it is necessary to isolate . for solution in the 

1 , le 

expression for I^, equation (3-133). To facilitate this, the terms not 

containing are combined in a parameter, such that 1^ becomes ^ 

- - 2it + + 2tan a^^[ln|sln - ln|cos(-^i^- 

- 2tan Oj^^Jlnlain -^^1 - ln|co8[-^^- a^g+i)|] (3-167) 


where Is given by 

le-1 ♦. ♦. 

l2tan aj[ln|slny-| - ln|cos (-j ®j)|] 

+ 2tan <»j^jlln|co8 “ ®i+l ^ I " ln|sln ^|]} 

+ I {2tan a [ln|sin -5-^1 - ln| 008(7=* - a ) | ] 
i-le+l ^ ^ ^ ^ 

+ 2tan [ln|cos (ji - ) I ” ln| sin ^| ] } (3-168) 

Similarly, Che expression for If,, equation (3-165), Is rewrlcten co isolate 
the unknowns K||, and u for solution* Thus, 


■c— ^ 




• f 
« • 
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where the constant terms are: 


a + 2N, 1/2 _ a + 2N, 1/2 _, ^ 

Cl - 4^s ' ^ i— t*"' ~) tan" -^] + 4Ky(^^,N^,b) 


(3-170) 

a + 2N, 1/2 _, a + 2N, 1/2 ♦ 

C 2 - 8i? - + 4(— ^ ) tan [ (— ^ -) tan -^] 


- 4Kpf6^,Nj,b) 

(3-171) 

- F 1 , il - 4v 

(3-172) 

*^ib " " ^ ^ 

(3-173) 


Using equations (3-123), (2-167), and (3-169), 1^ can be expressed as 
I4 - K^Cj + ♦ M C3 ^ - 0 (3-174) 

in which is defined, using the notation of equation (3-65), as 

J4 - aj bj Inde.le) ♦ Cj ln(le,l*<-l) ♦ ajlnlaln (3-175) 

where the coefficients are given by 

•j “ *^Ia ■*’ 


(3-176) 
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b, “ - 2 Can a, 

J le 


( 3 - 177 ) 


Cj - 2 can 


(3-178) 


Bj - 2 (can - can 


(3-179) 


An this point, the Integral constraint required for finite pressure 
gradients the trailing edge has been written with the unknowns isolated 
appropriately for solution. It remains to rewrite the previously given 
constraint equations of the Eppler method to include w as an unknown 
paraiseter. To aid in this, the quantity is defined by 


2w 

'wcl “ 

P ^ 

w 


(3-180) 


and Che equacion resulcing fron Che flrsc closure conacralnc, equacion (3-67), 
la wrlccen aa 


Si «cl ♦ 'Si «cl ^ ** 

P P 


(3-181) 


S«.l ♦ J.-O 


(3-182) 


't) 


where is defined as before, equaCions (3-69)-(3-73), excepC Chat Che 
expression given for a^, equacion (3-70), no longer includes Che Cere 
involving u. Siailarlyi Che equacio>. resulting froe Che second closure 
conscrainc, equacion (3-76), bscosws 
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where 


Wygi • J_ 

P 


(3-183) 


and Che u term is no longer Included In Che expression for equation 
(3-78)» which is used in defining J_, equation (3-77). Finally, equation 
(3-82), obCained using the Hatching condition at the trailing edge, is 
rewritten as 


In W^(0) + Kjj In W^(2e) - u In W^i2ir) +3^.-0 


(3-184) 


and, as before, the u tern is eliminated fron the definition of a^, equation 
(3-84), used in the expression for J^, equation (3-82). 

Equations (3-174), (3-181), (3-182), and (3-184) can be used to solve the 
four unspecified paraneters of the design problem. In matrix notation, this 
system of equations is represented by 


«=l 

W T 

cl 

p 

W , 
wcl 

1 

«sl 

W T 

sl 
. p 

W T 
V8l 
^ 1 

-lnW.(O) 

InW (28) 

-InWj 

8 

8 

w 

*=1 

^2 

^3 


(3-185) 

Expanding this system results in a transcsnden:.! equation for 4^^ which is 
! the form 


a ♦ b ln(le,le) ♦ c ln(le,levl) ♦ 


7 « 


e 



0 


( 3 - 186 ) 


As in the upnodlfied scheiae, this equation is solved by Regula falsi, after 
which, Jg, J,, and Jj are calculated and K^, K^, and u obtained by back 
substitution. At this point, it Is still necessary to satisfy equation 
(3-116) in order for the resulting airfoil to have a finite complex trailing 
edge pressure gradient. In order to achieve this, the computational procedure 
allows the designer to choose one of several possible Iteration schemes which 
vailes one of the parasMters, K, K, P, or until the condition of equation 
(3-116) is met. 


Calculation of Trailing Edge Pressure Gradients in the Modified Eppler Method 
Once an airfoil has been designed having bounded trailing edge pressure 
gradients, it is of interest to know the values of those gradients. As Zf" is 
obtained Ir ,*c?.-ction with the decign process from equation (3-98), then, in 
reference o :*on (3-28), the calculation of further requires 
that be determined. While z^^'' is made finite by requiring continuit y of 

1 V 

P**(9) at the trailing edge, the actual calculation of the value of by 

further differentiating the series representations of P(4) and <)(•»), in a 
manner similar to that used In obtaining ax'", would necessitate that P'O) 
be continuous everywhere and that P**(6) be piecewise contimious. As this 
would place severe restrictiona on the velocity distributions aliowable simply 
to facilitate the calcularlon of Rx> alcamaclve method was developed to 
approxlamte the value of The schema employed centers on expanding the 

tranaformstlon, s(c), in a Taylor series about the trailing edge to obtain 
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2(0 - 2^ + z^’(oc^) + It + 3T 

+ ;^ 2^'^ (C-Cj)^ + ... (3-187) 

In which, as determined previously, 

2 ^’ = 0 ( 3 - 188 ) 

* 1 (3-189) 

and, as just noted, z^" is given by equation (3-98). In addition, because the 
airfoil satisfies the conditions for a finite complex trailing edge pressure 
gradient, the value of zj" can be used in equation (3-27) to solve 
for By substituting a coordinate from near the trailing edge on the 

airfoil a.id the corresponding coordinate on the unit circle into equation 

1 V 

(3-187), is possible to approximate using the Taylor series 

reoresentation truncated to fourth-order terms. 

In ictual applications, by comparison with exact values obtained using 
airfoils from the von Mises family, it was found that the most reliable values 

i V 

of 2 ^ from Che approximate method were obtained when the upper surface 

coordinate nearest the trailing edge was substituted into the truncated series 

and this result averaged with one obtained using a similar point on the lower 

iv 

surface. A value of resulting from a obtained in this manner, 

however, must be viewed with some caution. The pr(.«ary difficulty is that the 
calculation is very sensitive to the detailed geometry of the trailing edge 
region. For example, while it odght be expected Chat Che shape of Che 
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trailing edge should not be f.ltered significantly by increasing the number of 
points used to generate the airioil, it can happen that when the resolution is 
improved by increasing the number of points, what appeared to be a reasonable 
trailing edge shape is actually overlapped. While this geometrical alteration 
might be so s' ight that it is not evident in plotted results, the influence on 
Rj is very significant. Although better accuracy can be achieved by greatly 
increasing the number of points used to generate an airfoil, this is entirely 
un 'cessary with regard to the basic design problem. Thus, in light of the 
large number of computer runs dictated by the iterative nature of design, and 
because the need for higher accuracy has yet to be established, the 
maintenance of low-cost and minimal run time were considered more important 
than obtaining Rf^ to higher accuracy. 

Finally, it should be noted that once the non-dimensional values of the 
trailing edge pressure gradients have been determined, they are resolved into 
components and denoted as Cpg<j^ and in the output from the modified design 

code. The streamline flowing from the trailing edge is considered to be 
directed along the bisector of the trailing edge closure angle and Cp^<p is the 
non-dimensional pressure gradient with respect to that direction. The 
non-dimensional pressure gradient normal to the trailing edge streamline is 
given by 

Influence of the Conditions for Finite Trailing Edge Pressure 
Gradients on Airfoils Designed Using the Eppler Method 

In order to better understand the Inpact of the finite complex trailing 
edge pressure gradient conditions on the design of airfoils using the Eppler 
and Somers code [25], It Is instructive 'co compar«> results obtained using the 
original code with examples generated using a version of the code which 
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incorporates the additional constraints. The airfoils to be considered in 
this comparison are only to aid in the understanding of the influence of the 
additional conditions on the designs generated by the code and are not 
necessarily intended as viable design possibilities. 

The airfoil and velocity distribution shown in Figure 25 were obtained 
using the unmodified Eppler and Somers code. In this extreme case, none of 
the iteration schemes for achieving a particular trailing edge angle are 
implemented and the trailing edge geometry resulting from the specified input 
produces very steep velocity gradients in the trailing edge region, as well as 
a very low trailing edge velocity ratio. Using the same input design 
parameters in the modified version of the code, the value of is determined 
by the method such that the integral constraint required for finite trailing 
edge pressure gradients, equation (3-II3), is satisfied. As seen in Figure 
26, imposing this constraint causes the flow in the recovery region and in the 
vicinity of the trailing edge to be modified considerably. In particular, 
note that the extent of the steepened gradients due to the closure 
contribution has been lessened and that the trailing edge velocity ratio 
increased significantly. This airfoil does not satisfy all of the conditions 
for finite trailing pressure gradients, however, in that the requirement of 
equation (3*-li6) has not been met. Next, the unmodified code is used with the 
same input as before, except now the value of K Is iterated to achieve ^3 ” 



iterated value of K has produced an airfoil very dissimilar to that obtained 
iisLng the value initially specified. Last in this series of comparisons. 
Figure 28 presents the airfoil obtained by Iterating the value of K from that 
used for the design of Figure 26 In order to satisfy equation (3-116). The 


global Influence, in particular the reduction in aft loading, which results 
from imposing the conditions for finite trailing edge pressure gradients 
should be noted. 

It must be emphasized that the class of airfoils having finite trailing 
edge pressure gradients is a subset of the family of airfoils designable using 
the Eppler code. Thus, the utility of the modified method is to facilitate 
the determination of the appropriate parameters which allow the conditions for 
finite trailing edge pressure gradients to be satisfied. If after these 
values have been found they are input into the unmodified! code, then the 
resulting airfoil would be the same as that generated by the modified scheme. 

As the airfoils considered in the preceding examples are somewhat 
extreme. Figures 29 and 30 provide a comparison of designs which are more 
reasonable. The airfoil shown in Figure 29 is obtained using the unaltered 
code while that of Figure 30 is a result of the modified version. As before, 
the reduction in aft loading and the increase in trailing edge velocity ratio 
occurring in the case of the airfoil generated with the additional constraints 
in force should be noted. 

Further appreciation of the behavior ot the modified code is obtained by 
comparing the differences In the manner that the modified and unmodified codes 
are used to design symmetrical airfoils. A symmetrical profile using the 
original Eppler and Somers code is obtained by setting corresponding upper and 
lower surface design parameters equal to one another. If iteration to a 
particular trailing edge closure angle is desired, a mode is chosen In which 
the selected upper and lower surface interating parameters are incremented In 
a manner that maintains the equality. Because the modified code solves for 
the value of u that allows the additional Integral constraint to be satisfied, 
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:tJk 

! 



its value cannot be specified equal to u. A symraetrial airfoil will result, 
however, if the other corresponding upper and lower surface inputs are set 
I equal, and the iteration mode selected to satisfy equation ((3-116) iterates 

on u- Fulfilling the requirement that w"(0) = w'*(2it) in this way forces u to 
j be equal to u. As demonstrated, for example, by the result shown in Figure 

r\ 31, this procedure yields a symmetrical airfoil having finite trailing edge 

* pressure gradients* 

1 

!. 

I 
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CHAPTER IV 


DESIGN EXAMPLES AND APPl. .ATIONS 


To explore the characteristics and capabilities of the modified version 
of the Eppler code, the usage of which is described in Appendix C, the 
velocity distributions of several airfoils appearing in the literature were 
adjusted as necessary and used to generate comparative airfoils having finite 
trailing edge pressure gradients. The first airfoil to be considered in this 
manner is the design of Strand presented in Reference [34]. The inverse 
method developed by Strand, used to generate this airfoil, is a development of 
Arlinger's procedure [10] which, in turn, grew out of that of Lighthlll [9], 

In the formulation of this procedure, the constraints imposed by the Inverse 
problem on the velocity distribution are satisfied by making adjustments to 
the portion of th ■ desireu velocity distribution which occurs on the lower 
surface of the airfoil. The differences between the desired distribution and 
that achieved are minimized by making the required adjustments as small, in a 
least~squares sense, as possible. The design point potential flow velocity 
distribution for this airfoil, shown in Figure 32, is calculated using the 
coordinates given by Strand [34] in the panel-method analysis procedure of the 
Eppler and Somers code [25]. At this angle of attack, the velocity 
distribution is intended to have a constant velocity rooftop followed by the 
appropriate Stratford recovery. The non-smooth appearance of the points on 
this calculated velocity distribution is largely a result of having an 
insufficient number of coordinates to describe the airfoil. Also shown in 
Figure 32 are the airfoil and velocity distribution obtained when the airfoil 
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is modified to have finite trailing edge pressure gradients. The two 
distributions are seen to be very similar except near the trailing edge. In 


fact, it was found that these oxfferences could be largely eliminated by 
splitting in additional coordinates. A comparison of the original Strand 
airfoil with the modified one is not presented as the geometrical differences 
between them are almost imperceptible. In comparing the analysis results for 
the two airfoils at the design angle of attack and Reynolds number of 3x10^, a 
lift coefficient of 1.32 is calculated for both airfoils. With natural 
transition determined by the program, the lift-to-drag ratio of the original 
airfoil at the design point is found to be 197, while that of the modified 
profile is slightly better at 207. 

Another example of a airfoil redesigned such that the trailing edge 
pressure gradients are bounded is that shown in Figure 33. The parent airfoil 
in this case is one developed by Llebeck and given the designation L1004 in 
Reference [12]. It is intended that this airfoil have a fully turbulent 
rooftop at the design Reynolds number of 3x10^. Although details of the 
velocity distribution and airfoil coordinates are unavailable, points taken 
from the velocity distribution presented in Reference [12] are noted in Figure 
33 for comparison with the distribution obtained for the modified version 
having finite trailing edge pressure gradients. While the design lift 
coefficients of 1.31, as calculated using the code of Eppler and Somers [25], 
is slightly less than the value given of 1.35 given for the L1004 by Llebeck 
[12], the calculated lift*to*drag ratio of 184 for the modified airfoil is a 
slight improvement over the value of 181 for the L1004. 

Another example, from the same family as the preceding airfoil, is the 
Llebeck L1003, de^-lgned to have a fully laminar rooftop at a Reynolds number 
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of 2x10^. As evident in Figure 34, the redesigned velocity distribution is 
close to that plotted in Reference [12]. In this case, however, even after 
fixing transition at the appropriate location, the Eppler and Somers analysis 
indicates that the entire recovery region is separated. As noted by Liebeck 
[34], obtaining an unseparated recovery for this section is extremely 
dependent upon having particular flow conditions at the beginning of the 
pressure rise. Thus, It is likely that the separation problems result because 
of small differences in the flow conditions between the two airfoils at the 
initiation of recovery. 

In addition to resulting in abrupt stalling behavior, the use of a 
Stratford pressure recovery for practical airfoils is often criticized in 
that, with regard to the normal flowfield variations that occur in realistic 
applications, the attainment of the precise flow conditions required at the 
beginning of the pressure rise cannot be assured. Consequently, the 
consideration of designs Incorporating Stratford distributions is often of an 
academic nature, as is the case of studies directed toward exploring ultimate 
possibilities in airfoil performance. Thus, an example of a more practical 
candidate for adaptation to an airfoil having finite trailing edge pressure 
gradients is provided by the well-proven and documented Wortmann FX 67-K-150, 
Reference [36]. The actual airfoil considered here, shown in Figure 35, is 
defined by the aerodynamically smoothed coordinates given by Somers [37]. 
Although this section is optimized for use with flaps, only the configuration 
having a neutral flap setting will be treated. The lack of steep gradients 
near the trailing edge in the velocity distribution examples considered up to 
now has permitted the conditions for a finite complex trailing edge pressure 
gradient to be applied to the existing designs with only minor alterations 



necessary* Hence, the goal has been to obtain an airfoil with bounded 
trailing edge pressure gradients that has a velocity distribution as close as 
possible to that of the parent section* The Wortmann section, however, is 
unlike those already considered in that, due to the significant differential 
between velocities on the upper and lower surfaces near the trailing edge, a 
steep adverse pressure gradient is present over the aft portion of the 
velocity distribution* Consequently, the concern in this case is to obtain an 
airfoil with finite trailing edge pressure gradients which, although having a 
somewhat altered velocity distribution, embodies the same design philosophy 
and achieves comparable performance* With this in mind, consider the result 
shown in Figure 37 and note that the use of the modified code has changed the 
velocity distribution such that the upper surface aft loading present on the 
Wortmann section, causing the steep adverse pressure gradient near the 
trailing edge, has been eliminated* It should be pointed out that there is a 
drag penalty associated with che steep lower surface favorable pressure 
gradient that Is a result of the closure contribution on the newly designed 
airfoil* While it was found that the drag could be reduced considerably by 
beginning the lower surface closure contribution sooner, and thereby lessen 
the gradient, the distribution shown was retained as it is more like that of 
the Wortmann section* In addition, in order to control separation problems 
that were introduced by the velocity distribution changes at the trailing 
edge, some modifications were made to the shape of the upper surface recovery 
distribution* Although the elimination of the upper surface aft loading 
results In a loss of lift as calculated by potential flow methods, this is 
more than offset by the increased value of Vj^/U which allows the lower 
surfact^ to carry a greater amount of aft loading* A comparison of the overall 
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aerodynamic performance of the two airfoils Is provided by the viscous 
analysis results obtained using the Eppler and Somers code and presented in 
Figures 36 and 38. The finite trailing edge pressure gradients airfoil has 
been designed such that the maximum lift-to-drag ratio occurs near the design 
angle of attack dictated by equation (3*111). Both airfoils exhibit best 
llft*to*drag ratios at a lift coefficient near unity. A more detailed 
comparison reveals that the drag polars of the new profile are roughly 
equivalent to those of the original section over roost of the usable 
performance range; however^ the performance of the section generated with the 
modified code is extended considerably in the direction of higher lift 
coefficients. 

While the imaginary part of the condition necessary for an airfoil to 
have finite trailing edge pressure gradients can only be satisfied at a single 
angle of attack, as equation (3*27) reveals, the results of the Eppler and 
Somers code viscous analysis of such airfoils indicates nothing particularly 
special about the aerodynamic characteristics at that angle of attack. This 
fact, however, should be to some degree expected in that the analysis makes 
use of conventional boundary layer theory in which normal pressure gradients 
through the boundary layer, as well as all wake influenceb, ^>*e assumed to be 
unimportant. As has been discussed, this assumption breaks down near the 
trailing edge where the inviscid pressure gradients are generally unbounded. 

It should be no surprise, then, that the results of a calculation based on 
conventional boundary layer theory do not indicate any characteristics 
attributable to the presence of finite trailing edge pressure gradients* 

Thus, a thorough evaluation of the effect of such factors on airfoil 
performance would require a fairly extensive investigation that makes use v>f a 
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theoretical model having a more detailed description of the flow in the 
vicinity of the trailing edge. To demonstrate what such a model might 
indicate, consider the application of the GRUMFOIL code [21] to the two 
airfoils just presented. The results of this malysls for the FX 67-K-150 are 
given in Figures 39-41, and those for the corresponding finite trailing edge 
pressure gradient airfoil in Figures 42-44. The analysis of the latter 
section was performed at its design angle of attack and that of the Wortmann 
at an angle of attack which resulted in the lift coefficients matching. The 
Reynolds number used was 2x10^ and the Mach number was set to zero. The 
aerodynamic characteristics calculated using GRUMFOIL are somewhat different 
from those obtained with the Eppler and Somers code. In general, the lift 
coefficients calculated by GRUMFOIL are slightly greater than those of the 
Eppler and Somers code, while the drag coefficients, even though transition 
predictions agree fairly well, are notably less. 

The fully viscous GRUMFOIL pressure distributions for the two airfoils 
are given in Figure 39 and 42. In considering the viscous pressure 
distribution for the Wortmann section. Figure 39, the pressure spike near the 
leading edge which is not present in the Eppler and Soibors potential flow 
results warrants explanation. Based on distributions obtained at lower angles 
of attack, it was concluded that the peak is due to a lack of smoothness in 
one of the coordinates rather chan from the angle of attack under 
consideration being too large. It has been found that only a very small 
inconsistency in the given coordinates can be responsible for such a result. 

In further considering Che viscous pressure distribution of Che Uortiaann 
design, it should be observed that Che steep upper surface gradients near the 
trailing edge present ir the potential flow results are largely eliml- ^ed by 
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the smoothing action of viscous influences. As seen in Figure 42 for the 
airfoil having finite trailing edge pressure gr^dients, on the other hand, the 
inviscid calculations are little impacted by the inclusion of viscid-inviscid 
iterations. 

Further differences in the trailing edge region flow behavior are well 
demonstrated by comparing the boundary li^yer characteristics of the two 
airfoils. In the case of the Wortmann profile, Figures 40 and 41 readily 
demonstrate the singular behavior at the trailing edge of the displacement 
thickness, form factor, and the equivalent surface source velocity. It should 
be noted that these results conform very well to those found experimentally, 
such as In the work of Preston, et. al., (29]-(31). In remarkable contrast, 
as seen in Figures 43 and 44, the slope discontinuities are eliminated for the 
airfoil with bounded trailing edge pressure gradients. The ramifications of 
these results are significant. In addition to any performance benefits 
arising from smooth flow off the airfoil and into the wake, the application of 
the condition for finite trailing edge pressure gradients has produced a class 
of airfoils for which the strong viscld-lnvlocld interactions, bcyr.id that of 
the displacement thickness, an be neglected. That Is, convt' • ' • ial boundary 
layer theory remains valid in the region of the trailing eo ''•? -’uch 
airfoils and is sufficient for the prediction of their aecoO./ -1.. 
characteristics. Furthermore, as the Influences die to vlscosi'v / 
minimised, the results calculated using potential flow design shoul 1 

be more reliable than those generally obtained. 

Because the imaginary part of the condition required for achieving finite 
trailing edge presaure gradients can only be satisfied at a single angle of 
attack, equation (3-27), it la of Interest to examine the ieportance of this 
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limitation by consUering the off^deslgn boundary layer behavior as calculated 
using GRUMFOIL. First, however, it should be noted that becaus. of 
differences in the zero-lift angle of attack prediction, it is likely that the 
correspondence between the angles of attack calculated by the Eppler and 
Somers code and those used in GRUMFOIL is not exact* Thus, in all 
probability, the case already presented represents a slightly off-design 
situation. In any event, to further consider the flowfield behavior off- 
design, Figures 45 to 47 sumouirize the GRUMFOIL output for the Wortmann based 
finite trailing edge pressure gradients airfoil at an angle of attack of 
approximately four degrees less than the design value. From these results, it 
is apparent that strong singular boundary layer characteristics at the 
trailing edge do not dramatically appear when the airfoil is operated at 
conditions other than those of the design point. It is evident f'^om Figure 
45, however, that the streamwise change of pressure along the wake centerline 
in the vicinity of the trailing edge, and consequently its effect on the 
Inviscid flow, has increased over chan seen in Figure 42. In addition, 
although certainly not discontinuous, the slopes in the immediate vicinity of 
the trailing edge on the boundary layer property distributions do appear 
slightly steeper than those present at the design angle of attack. Uthough 
additional verification is warranted, on the basis of these off-design 
GRUMFOIL results, it can be concluded that if any aerodynamic benefits are 
realized by the pvesence of finite trailing edge pressure gradients, Chen 
these henefi*^a are not liodted to the design angle of attack but are present 
to some extent er an operational range of angles. Thus, in addition t^ 
being of academic interest, this allows airfoils designed with finite trailing 
edge prsssurs gradients to marie considsrstion for practical application. 


CHAPTER V 


C0NCLUDK4G REMARKS 


The potential flow sol iClon for any airfoil having non-'zero trailing edge 
loading is characterized by the presence >.>f unbounded pressure gradients at 
‘'he trailing edge. Although in a real fluid the pressure gradients are 
somewhat softened by viscous affects, those in the trailing edge region do 
* ideed become extremely steep resulting in, among other things, the slop*, of 
the displacement thick. »» distribution being disconrir.uoun it the t. railing 
edge. Considering the near critical nature of eaay of the velocity 
di‘'.tributions prescribed for maximum lift or minimum drag, the en‘*ounter of 
such a disturbance could be sufficient to cause severe upstream sepa.-ation 
problems. Thus, the goal of separation free flew should benefit by the 
removal of this disturbance to allov the fluid on the airfoil to flow into the 
wake as smoothly as pcssible. Of additional concevn in this regard, the 
presence of strong adverse pressure gradients in the vicinity of the trailing 
e<ge, as seen in many maximum performance design efforts, nay result in 
upstream separation problems. Thus, to help ensure that the hlgii performance 
levels promised by potential flow methods arc realized in practice, a 
procedure haa been developed to design airfoils for which the trailing edge 
pressure gradients are finite and the flows on the upper and lower surfaces 
approach the trailing edge ree of strong adverse pressure gradients. 

The ability to specifically configure the trailing edge region of an 
lirfoil to achieve finite pressure gradients has been made possible by the 
unique capability of the Eppler method which allows different segments of an 
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airfoil to be designed for different angles of attack. As the removal of the 
trailing edge pressure gradient singularities requires that no load be carried 
by the trailing edge, the method is able to adapt the aft portions of the 
airfoil such that this no-load requirement is met. To some extent, the 
resulting airfoil can be thought of as one in which the trailing edge region 
behaves locally like a flat plate at zero angle of attack although, in the 
flat p’ate case, gradients are not only finite but zero* By eliminating 
the unbounded trailing edge pressure gradients, it is possible to specify a 
velocity distribution on an airfoil which pushes boundary layer performance to 
its critical limits as is the case, for example, in specifying a Stratford 
recovery, or in choosing a distribution for which the trailing edge velocity 
ratio is maximized* 

In the formulation of conventional boundary layer theory, normal pressure 
gradients through the boundary layer are ignored and only the influence of the 
displacement thickness on the invlscid results is considered* Thus, because 
of the unbounded pressure gradients that generally occur at the trailing edge 
in the potential flow solution, conventional boundary layer theory is invalid 
in the vicinity of the trailing edge* In regard to this limitation, Melnik, 
et* al* [33] demonstrated that the potential flow solution singularities give 
rise to additional vlscid**inviscid interactions, each having an effect as 
Important as that of the displacement thickness* Thus, by allowing for the 
influences caused by the normal pressure gradients in the trailing edge 
region, wake thickness, and wake curvature, Melnik and his coworkers developed 
a self-consistent boundary layer theory able to account for the strong viscous 
interactions due to the singularities in the inviscid flow solution* Although 
the formulation is distinctly different, the removal of the trailing edge 
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singularities can be considered an alternative approach to the same problem. 

In this light, airfoils having finite trailing edge pressure gradients 
represent a class for which the strong viscid-inviscid interactions in the 
trailing edge region have been minimized. Consequently, conventional boundary 
layer theory is sufficient for the viscous analysis of such airfoils. 
Furthermore, because the corrections necessary to the inviscid solution due to 
viscous effects are minimal, potential flow design methods are likely to yield 
more reliable results than they otherwise would. 

Considering the nature of the flow behavior in the region of the trailing 
edge, airfoils designed to have finite trailing edge pressure gradients may be 
ideally suited to aid in the development and calibration of improved 
aerodynamic prediction methods for airfoils. For example, in the theoretical 
formulation used in the GRUMFOIL code [33], the local trailing edge region is 
modeled as unseparated flow over a flat plate at angle of attack. Thus, the 
class of airfoils having bounded pressure gradients at the trailing edge are 
much more consistent with this model than is generally the case. Such 
airfoils should, thereiore, provide useful development tools and calibration 
cases. In a similar application, because the rapid growth of the displacement 
thickness at the trailing edge that generally occurs leads to numerical 
divergence problems, the development of viscous analysis methods in which 
potential flow-boundary layer iteration is employed should benefit from the 
well-behaved growth in displacement thickness at the trailing edge on airfoils 
having finite trailing edge pressure gradients. 

Finally, if imposing the requirement for finite trailing edge pressure 
gradients does Indeed minimize the viscous interactions and allow potential 
flow predictions to be more fully realized, then this situation clearly 
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suggests that an improvement in airfoil performance is possible. In the 
examples that were considered » the design effort was directed at matching the 
characteristics of previously defined velocity distributions. Consequently, 
it remains to explore the potential of exploiting the use of finite trailing 
edge pressure gradients to enhance airfoil aerodynamics. Encouragement that 
gains might be made, however, is provided by the results of the GRUMFOIL 
analysis from which, for example, the employment of the finite trailing edge 
pressure gradients condition yields an airfoil having a thinner displacement 
thickness and wake than otherwise occurs. If such performance benefits are 
Indeed found to exist then, as the GRUMFOIL results indicate that reasonable 
off-design capability Is present, airfoils having finite trailing edge 
pressure gradients become candidates for practical application. 
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APPENDIX A 

LIMITING COEFFICIENT VALUES OF THE GENERAL TRANSFORMATION 
FOR MAPPING A CIRCLE TO M AIRFOIL 

Consider Che general transformation which maps a circle centered at the 
origin of the C-pIane into an airfoil in the z^plane as given by 


(t) 


1 - 

n-l C 


(A-1) 


where c^^ » + ib^. In this transformation, depicted in Figure 48, it is 
assumed that jc| > r and r > 1. G^ven that the origin of Che z-plane lies 
within the profile, the area enclosed by the boundary of the airfoil is 
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and, from the Cauchy-Rlemann equations. 
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Substitution of equation (A~3) Into (A-2) yields 
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the integral in equation (A-A) yields 
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Thus, equation (A-4) becomes 
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which can be rearranged to give 


A 

K 


« • n c 

,2 - J ' 


n»l r 


,2n 


(A-5) 


(A-6) 


d4i 


(A-7) 


(A-8) 


Since the area of the airfoil cannot be negative, It follows that 


ORIGINAL PAGE 13 
OF POOR QUALITY 


98 


> I 


“ n c 


n»l r 


2n 


(A-9) 


As the expression is valid for all values of r between one and infinity, it 
must hold for r equal to one. Thus, 
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An immediate consequence of this result is that 
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To examine these results further, consider the case for which the 
equality in equation (A-ll) holds, i.e., when |cj| is unity. Observe from 
equation (A-IQ), that for this to be true all of the other transformation 
coefficients must be zero. Thus, for this case, the transformation becomes 
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where |c|| * 1* Writing the transformation coefficient as 
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where 0 < y < 2ir, the mapping funcclon becomes 
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For the case of mapping the unit circle, th^s gives 

- 2.*’' co.(*-r) (A-14) 

Hence, as shotm in Figure 49, if ( describes Che unlc circle chen its 
conformal Image, z, describes boch sides of a flat place orlenced to the 
poslClve real axis aC che angle y. Thus, |cj| can be equal Co unicy onl/ for 
funcclons mapping che unlc circle Co a flac place. Noce chaC wlchouC loss of 
generallcy, che cralllng edge of che profile generaced can be assumed co be 
locaCed on che real axis. For che example given, Chls resulcs In Che 
oriencacion of che flac place being along che real axis wlch y - 0* Thus, a^ 
is unlCy and all che oCher a^'s and b^'s muse be zero. This case is 
equlvalenc co chaC of che Joukoiwsky flaC place aC zero angle of aCCack. 

To examine Che relaclonshlp of che Cransf ormacion coefflclenCs co che 
maximum possible cralllng edge velocicy, VTE, consider che expression for che 
Cralllng edge veloclCy for an airfoil obCalned from a mapping of che unlc 
circle as given by 
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where, for che case of che cralllng edge fixed on che real axis. 
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Izx'l ” ^ n(n+l)a + ( I n(n+l)b (A-16) 

n-1 " n-1 " 

Clearly, the maxiimim value of VTE/U is obtained when |z^"| is minimized. This 
occurs when the second term on the right is zero as accomplished when all of 
the b^'s are zero. Although there are non-zero values of the b^^'s which can 
achieve the same result, there are none which can result in a higher value of 
VTE/U. Consequently, the symmetrical airfoil that results when the b^'s are 
/.ero has a trailing edge velocity ratio which is at least as great as any non- 
syimnetrlcal airfoil generated using the same set of a^'s. 
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APPENDIX B 


COEFFICIENTS OF THE INEQUALITY EXPRESSION FOR POSITIVE THICKNESS 


The expression which is developed to Insure chat physically realizable 
airfoils will result from a six term von Mises transformation, equation 
(2-58). is 

Af^ + + Cf^ + Df + E > 0 (B-l) 

the coefficients terms for this expression are given by 

A - 16r^(r-l)^ (B-2) 

B • 8r^[5(r-l)^ - + a 2 (r-l)^ - a^Cr-l)^ 

♦ *4(r-l) - aj] (B-3) 

C ■ 4r^[(l0 - 4aj - 2a2)(r-l)^ + (3a2 ♦ 3a^ •** ajr)(r-l)^ 

- (2aj ♦ 4a^ + 2a^r)(r-l) > (a^ 5aj ♦ 3ajr)] (B-4) 

D ■ 2r[(l0 - 6 aj - 6 a 2 - 3aj) (r-l)^ ♦ ( 3 a 2 6 aj * Ba^ 

♦ 2ajr ♦ 4a^r) (r-l) - (aj + 4a^ + lOaj 2a^r 

lOa^r ♦ (B-5) 
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APPENDIX C 


USAGE OF THE EPPLER CODE INCORPORATING THE CONDITIONS FOR FINITE TRAILING 
EDGE PRESSURE GRADIENTS AND LISTING OF PROGRAM MODIFICATIONS 


In asking modifications to the Eppler and Souiers code 125] in order to 
facilitate the design of airfoils having finite trailing edge pressure 
gradients, the effort was made to leave as much of the existing code and its 
data input as unchanged as possible. Thus, the discussion included in this 
appendix should be considered in conjunction with the code description and 
operating instructions presented in Reference (25). 

While the primary purpose of the modified code is the design of airfoils 
having finite trailing edge pressure gradients, it might also be of use in the 
design of airfoils in which the upper and lower velocity distributions merge 
smoothly at the trailing edge without the nearby presence of steep gradients. 
In either case, the integral constraint of equation O-llS) is satisfied. 
Because this condition eliminates the pressure loading at the trailing edge, 
the shape of the aft portion the airfoil is largely governed by the zero 
closure angle which results. Thus, although control remains over the extent 
that this zero closure angle is allowed to influence the overall shape of the 
rear of the profile, much of the ability to iterate to a desired closure angle 
that Is present in the original code is lost. 

With reference to Eppler and Boners (25), the input to the modified code 
differs fron that of the original as follows: 

1. On the TRAl card, although assuaad internally by the modified code, 

o« should be set equal to a, as required by equation (3>132). 

P ‘ 
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2. The value of u is no longer specified by Che word on Che TRA2 
cardt buc decermined by Che program such chac Che fourch InCegral 
conscrainc is sacisfied. In ics place, however, Che quancicy IWPPM 
Is deslgnaced Co selecC Che mode of iceraclon used Co achieve Che 
requiremenc ChaC P”(0) • P”(2n), equacion (3-116). The iCeraCion 
mode posslbillcies are as follows: 

IWPPM *0 - No iCeracion is performed and P"(0) will, in 

general, noc be equal Co 
IWPPM - 1 - K is replaced by K ♦ &K 

IWPPM - 2 - K is replaced by K + AK 

IWPPM "3 - w is replaced by u + An 

T* —* 

IWPPM *4 - A is replaced by X + AA, unless chaC resulc is 

calculaced co be less Chan zero or greacer chan X. 
In chac case, Che program swlCches Co IWPRl * 1. 

3. Also on Che TRA2 card, Fg - RNSis> which dccerulnes Che 
incerprecacion of F 9 and P^q, ausC be sec co zero. Consequencly, 

.1 P 9 is always incerpreced as K. 

4. While scill acCive In Che code, ic should be noced chac Che 
specif icaclon of ITMOD equal co S or 6 , P|| on Che TRA2 card, will 
generally noc resulc in convergence Co Che specified value of K^. 
This is because in chase modes, K is iCeraCed and Che calculaced 
ICeracion incremtne, AK, is superseded by Che code decermlnacion 
of h which, in Curn, alters Che value of K* 


A sample input set for Che modified code is prcsenced in Figure SO. The 
resulc of chis input is the eirfoil having finite trailing edge pressure 
gradients shown in Figure 37. 
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Ic should be noted Chat chat It is not Intended that Che modified code be 
run with ITMOD * 0 and IWPPM * 0 simultaneously. The design of airfoils 
having finite trailing edge pressure gradients is carried out with ITuJD ^ 0 
and IWPPM * 0, The design of airfoils In which it is only desired that 
velocity distribution be free of steep adveise gradients In the vicinity of 
the trailing edge can be accomplished with ITK- <j * 0 and IWPPM ■ 0. In 
designing airfoils having finite pressure gradients, it Is sometimes 
advantageous to begin the process with one of the ITMOD isodes and switching to 
one of the IWPPM options when close to the desired velocity distribution. For 
example, a symowtrlcal airfoil can be obtained by first setting ITMOD ■ 6 or 
9, IWPPM • 0, and specifying Che additional inputs as described * * Cppler and 
Somers (23] such that the upper and lower surface velocity distribution 
specification quantities arc equivalent. In the case of using the modified 
code, however, the program will solve for a value of v which is different from 
that specified fer w. Now, inputting the results from this run and changing 
to ITMOD • 0 and IWPPM * 3, p will be iterated until it agrees with p and a 
syasMtrlcal profile having finite trailing edge pressure gradients will 
result. 

A listing of Che modifications made to the Epplcr and Somers code (23] 
will follow. Only the main program and modified subroutines, in addition tc 
several newly added subroutines, will be presented. The reader is again 
referred to Rcfcrenea [25] for a listing of Che original code. 
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PROGRAM EPPLER3 (INPUT, OUTPUT, TAPE4,TAPE5=INPUT,TAPE6=OOTPUT) 
DIMENSION XF(121) ,TF(121),BETAF(121) 

DIMENSION AM(7,7),AV(7) 

DIMENSION V( 14) ,MARKEN(20) ,ALCA(14) ,CAE(2) 

DIMENSION RE(5) ,MA(5) ,MU(5) ,T(42) 

DIMENSION TM(5) ,ALS(5) ,RER(5) ,MUR(5) 

DIMENSION TST(5) ,BANT(5) ,CW(5,2, 14) ,SU(5,2, 1~>) ,SA(5,2, 14) 

COMMON P1(121),P(121. ,XP(121),YP(121) ,PUFF(14) ,AGAM( 14) ,X( 121) , 
1Y(121) ,DS( 122) ,VF( 121) ,'ARG( 121) ,ANI(28) ,ALFA(29) ,IZZ,KFU,NQ,NUPRO, 
2JAB,JST,CM,ETA,ABFA,PI,BOGEN,DARG,PURES(13) ,FUW(60,7) ,RS(60) 

COMMON XTF,SMA,XFL(10) ,GAMMA(121,2) ,AMAT(120, 120) 

COMMON 'GRZK/CDK,AA(7) ,BB(7) 

COMMON/ PRAL/DLT , DLTU , ALN ; ALV ( 1 4 ) , NAL , ITP , NAMP ( 1 2 ) ,CML ( 1 4 ) ,CRL ( 1 4 ) 

1 ,CPV(2) ,ALn(4,2) ,DARF,ITIT1 ,ITIT2 
COMMON/ PLIM/MPL ,MGC ,XZEH , YZEH ,MSPLI 
COMMON/EA/ILES, IDRU ,ISTA,NNESE 
COMMON/TRIT/DLV,SUMP,XTRI(4) ,NU,ND 
COMMON /LINING/BR0KL( 12) ,NLINE(5) ,NPAR(5) ,JNEW 
EQUIVALENCE(XF( 1 ) ,FUW( 1 , 1 ) ) , (YF( 1 ) ,FUW(3, 3) ) , (BETAF( 1 ) ,FUW(5,5) ) 
EQUIVALENCE (CW( 1 , 1 , 1) ,P1 ( 1 ) ) , (SA( 1 , 1 , 1 ) ,P(20) ) , 

1 (SU(1,1,1),XP(39)) 

DATA ILES,IDRU,ISTA,NNESE/5,6,4,1H1/ 

DATA MARKEN/ 4HTRA 1 , 4HTRA2 , 4KALFA , 4HAGAM , 4HABSZ , 4HSTRK , 4HENDE , 
14HDIAG,4HR£ ,4HSTRD,4HFLZW,4HPLWA,4HPLW ,4HTRF ,4HAPPR,4HCDCL, 
24HPAN ,4HFXPR,4HFLAP,4HPUXY/ 

DATA CPV/9HVELOCITY ,9HPRESSURfc /,ALTX/4H ZER , 4H0-LI , 4HFT L, 

•4HINE ,4h CH0,4HRD L,4HINE ,4H /,KBLT/1H / 

DATA MGC , ISTIFT,MXZ ,CDK/0, 1 ,- 1 , . 0 1 / 

DATA ZAEH,DICHTE/13.6E-6,.12533/ 

MPL=0 

PI = 3.141592654 
BOGEN = 0.0174532925199 
ABFA = 1.0 
AGAM(2)=1. 

AGAM(3)=1. 

AGAM(6)=1. 

AGAM(8)=0. 

AGAM(10)s1. 

9 MTRsO 

1 1 READ(ILES,2)MARKE,NUPA,NUPE,NUPI,NUPU,PUFF 

2 F0RMAT(A4,3I1,I3,14F5.2) 

DO 12 1=1,20 

IF (MARRE.EQ.MARKEN(D) GO TO 13 

12 CONTINUE 

14 WRITE(IDRU,3) MARKE 

3 FORMAT (11H INCORRECT ,A4,5H CARD) 

GO TO 11 

C TRA1TRA2ALFAAGAMABSZSTRKENDEDIAGRE STROFLZWPLWAPLW TRF APPRCDCL 
C PAN FXPRFLAPPUXY 

13 GO 10(15,22,333,14,142,90,150,104,30,112,30,71,60,14,14,160, 
•170,180,190,106),! 
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THAI CARD 

15 NUPRO=NUPU+1000»MUPI 
IF(MTR.EQ.O)JSTsO 


1=0 

18 1 = 1+1 

ANRI=RUND(PUFF(I)*ABFA, 1000.) 
IF(ANRI.NE.O.)GO TO 20 
IF(JST.NB.O)GO TO 21 


JST=MTR+1 
20 MTRaMTR+1 
ANI(MTR)=ANRI 


1 = 1+1 

ALFA(MTR)=PUPF(I) 
IF(I.NE.14)GO TO 18 

21 JAB=MTR 
GOTO 11 

C TRA2 CARD 

22 DO 23 1=1,13 

23 PURES(I)sRUND(P0FF(I),1000.) 
MSPLIsO 

ITPsO 

IZZ=INT(PUFF(14)) 

CALL TRAPRO 
XDAsO. 

YDAsO. 


DEFLG=0. 

GO TO 9 
C RE CARD 

25 IF(PUFF(2).EQ.O.) GO TO 28 
DO 27 J=1,5 

RERX= PUFF(2»J) 

IF(RERX.EQ.O.)GO TO 26 

RE(J)=1.E5*RERX 

IPU = INT(P0FF(2«J-D) 

MA(J) = IPU/100 

MU(J) = IPU/10 - lOniA(J) 

27 JR = J 

26 DO 29 J=1,4 

29 XTRI(J)=PUFF(J+10)».01 

28 CALL GRP(NAL,RE,HU,JR,1STIFT) 

HSPLI=0 

JP=JR 
GOTO 11 
C FLZU CARD 

30 IF(NUPA.EQ.O) GO TO 41 
AGAM(6)«FL0AT(N0PE) 

AGAM(8)=PL0AT(NUPI) 

31 1F(I.EQ.9) 00 TO 25 
IF(PUFF(2).BQ.O.) GO TO 36 
GDF = PUFF(I) ■ 

VMAX a PUFF(2) 

IF(PUFF(3).ME.O.) DICHTE a . l«PUFF(3)/9.806 
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IF(PUFF(4).NE.0.)ZAEH = PUFF(4)»1 ,E-6 
IF(PUFF(5).EQ.O.) GO TO 50 
D = 0. 

DO 34 J = 1,5 
JZ = 2»J + 3 

IF(PUFF(JZ).EQ.O.) GO TO 36 
TM(J)=PUFF(JZ) 

ALS(J)=POFF(JZ+l) 

MUR(J) = NUPU 
34 JT = J 

36 IZT=NZPZ(2,6«NAU2) 

JP=JT 

WRITE(IDRU,37)IZT,NAMP, (ALTX(J,ITIT2) ,J=1 ,4) 

37 FORMAT (A1,36HAIRCRAFT-0RIENTED SUMMARY AIRFOIL ,12A1,3X, 
•31HANGLE OF ATTACK RELATIVE TO THE,4A4) 

I VMAX= INT ( VMAX* 3 . 6 ) + 1 
IZT=NZPZ(2,0) 

WRITE( IDRU , 38 ) IZT ,GDF , IVMAX ,DICHTE , ZAEH 

38 FORMAT (A1,6H W/S =,F6.2,8h KG/SQ.M,3X,7HV MAX =,I4,5H KM/H.3X, 
•5HRHO =,F5.3, 13H KG«S E2/M E4,3X,4HNU =,F10.8,7H SQ.M/S) 

IZT=NZPZ(2,0) 

WRITE(IDRU,40)IZT, (KBLT,TM( J) ,ALS( J) ,J=1 , JT) 

40 FORMAT (A1 ,5X,5(A1 ,4X,3HC =,F5.2,8h THETA =,F5.2)) 

41 VI s SQRT(2.»GDF/DICHTE) 

DO 48 1=1, NAL 

IVS = -I 

DO 46 J s 1 , JT 

IF(ALV(I)-ALS(J))42,44,42 

42 VALF = V1/SQRT(.11»ABS(ALV(I)-ALS(J))) 

IF(VALF - VMAX)46,46,44 

44 VALF = VMAX 

46 RER(J) = VALF»TM(J)/ZAEH 

48 CALL GRP(IVS,RER,MUR,JT,ISTIFT) 

MSPLIsO 

49 IF(D)72,50,72 

50 GO TO 11 
C PLW CARD 

60 IF(PUFF(1) >62,68,62 
62 DST s PUFF(1)».01 
GST = PUFF(2) 

DGF = PUFF(3) 

CWSF s PUFF(4)».001 
DO 66 J 3 1,5 
JZ s 2«J-».3 

IF(PUFF(JZ) >64,68,64 
64 TST(J) s PUFF(JZ) 

BANT(J> s PUFF(JZf1> 

MUR(J> = NUPU 
66 JT = J 
66 BF s 0. 

FST s 0. 

MF * 0 
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CALL DIA(X,Y,NQ,D) 

IFvDST.LT.O.) DST = D 
DO 70 J s 1,JT 
TM(J) = TST(J)»DST/D 
ALS(J) 3 0. 

BF s BF + BANT(J) 

70 FST s FST + BANT(J)»TST(J) 

FF = FST*DST/D 
OEW = GST + (FF-FST)»DGF 
GDF s GEW/FF 
GO TO 36 

C PLWA CARD 

71 NAN 3 NUPU 
CWSFUsCWSF 
FAU 3 FF 
GAU 3 GEW 
DO 88 NF 3 1,NAN 
FF 3 FF+PUFFd) 

GEW 3 GEW4>PUFF(2) 

CWSFsCWSF+ . 00 1 •PUFF ( 3 ) 

VI 3 SQRT(2.»GEW/(DICHTE«FF)) 

72 I2T3NZPZ(3,NAL+10) 

WRITE( IDRU , 74 ) IZT , NAMP 

74 FORMAT (A1 ,25HAIRCRAFT POLAR AIRFOIL ,12A1) 

CWS 3 CWSF/FF 
IZT3NZPZ(2,0) 

WRITE(IDRU,76)IZT 

76 FORMAT (A1,39H B(M) S(SQ.M) S«(SQ.M) W(KG) W*(KG) 
•3X , 3HT/C , 3X , 6H ( T/C ) • , 2X , 8HAP ( SQ .M ) , 2X , 3HCDP ) 
IZT3NZPZ(1,0) 

WRITE( IDRU , 7 8 ) IZT , BF , FF , FST ,GEW »GST ,D ,DST , CWSF ,CWS 
78 FORMAT (A1 ,F6.2,2F8.2,2F8.0,4F8.4) 

IZT3NZPZ(2,0) 

WRnE(lORU,80)IZT 

80 FORMAT (A 1,54H ALPHA CL CDP CDT V(KM/H) 
•/D) 

DO 84 I 3 1,NAL 
CA 3 0. 

CWP 3 0. 

DO 82 J 3 1,JT 
BATF3BANT ( J ) •TSTC J ) /FST 
CALL Vised J,C ANT, CWNT,CMDU) 

CA3CA4>CANT«BATF 
82 CHP3CWP4CWNT«BATF 

CHGES 3 CWP ♦ CWS ♦ 1.03»CA»CA»PF/(PI»BP»BF) 
IP(ABS(CA) .LT. .ODCAa.OI 
VKMH 3 3>6«V1/SQRT(ABS(CA)) 

VS 3 VKMH*CW0ES/(3.6*ABS(CA)) 

GLTZ 3 CA/eWGES 
IZT3NZPZ(1,0) 

84 WRITE(IDRU,86) IZT,ALV(I) ,CA,CWP,CWQES,VKHH,VS,QLTZ 
86 FORMAT (A1 ,F6.2,F8. 3,2F8.4,F8. 1,F9.3,F8.2) 


E 


1 




j 

} 
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(t) 


IF(NF.EQ.O) GO TO 50 
88 CONTINUE 


CWSF = CWSFU 
FFsFAU 


GEW = GAU 
GO TO 50 
C STRK CARD 

)0 IF(NUPU)94, 100,92 

92 NT = 0 

94 DO 98 J = 1,14 

IF(PUFF(J)) 96,100,96 
96 NT = NT + 1 
f8 T(NT) = PUFF(J)«10. 

IF(IABS(NUPU).GT. 14) GO TO 11 
no CALL STRDR(T,NT) 

IF(NUPI.NE.O) GO TO 11 
DO 102 I s 1,NT 

102 CALL STRAAK(T(I),RUA,YBL,MXZ,ISTIFT) 

GO TO 11 
C 01 AG CARD 

10<4 CALL DIAGR(ISTIFT,NUPU,NUPI) 

GO TO 11 
C PUXY CARD 
106 CALL PUDECK 
GO TO 11 
C STRD CARD 

112 IF(NUPU.NE.O) MXZ s NUPU 

IF(PUFF(1).NE.0.) YBL = 100.»P0FF(1) 

IF(PUFF(2).NE.O.; RUA = 100.«PUFF(2) 

GO TO 11 
C ABSZ CARD 

142 IF(NUPA.NE.O) AGAM(3)sFLOAT(NUPE) 
i:?(PUFF(2).NE.O.) ABFAsPUFF(2) 

GO TO 11 
C ENDE CARD 

150 IF(MGC.NE.O) CALL GCLOSE 
IF(MPL.NE.O)CALL FINISH 
STOP 

C CDCL CARD 

160 IF(NUPA.EQ.O)G0 TO 166 
BL1sPUFr'J)-.005 
BL2ari?F(2)+.005 
DO 162 Ka1,5 

?!LINE(K)aINT(BL1*(10.»*(K-3)))-10»INT(BL1»(10.«»(K-4))) 
le-* MPAR(K) aINT(BL2»(10.»»(K-3)))-10»INT(BL2»(10.»»(K-4))) 
DO 164 Ks1,12 
1 64 BR0KL(K)*PUFF(K+2) 

GO TO 11 

166 CALL CDCL ( NUPU, JP,ISTIFT) 

GO TO ;i 
C FXrS CARD 

180 :tp*nupu 


i 


I 


i * 
• I 




i 


I 


I 

) 

I 


I 

I 


t 

i 


I 


! 


I 

I 

I 



t 



original page 

OF POOR QUALITY 


111 


CALL PULES 
MSPLIsO 


C PAN CARD 

170 IF(MSPLI.EQ.O)CALL SPLITZ(X,Y,NQ,XP) 

IF(MUPA.NE.O) AGAMdO)sFLOAT(NUPE) 

IF(NUPA.EQ.9) GO TO 11 
DO 171 Is1,14 
IF(PUFF(I).EQ.O.)GO TO 172 
MEIG=INT(PUFF(D) 

MEI=MEIG/ 1 0 

KEl ±MEIG-10«ME1 

XSTXsABS(PUFF(I)-FLOAT(MEIG) ) 

CALL PADD(X,Y,XP>NQ,MEI,KEI,XSTX) 

172 CONTINUE 

DO 174 Is1,NQ 
XF(I)»X(I) 

YF(I)=Y(I) 

174 BETAF(I)=XP(I) 

DLTRsDLT 

DLTURsDLTU 

XDAsO. 

YDAsO. 

FLCHsO. 

DBFLQsO. 

NQRSsNQ 
176 NKRsNQ-1 

CALL PANEL(NKR|AMAT,GAtMA,CAE) 

DARG = ALN 
GO TO 11 
C FLAP CARD 

190 CHORD s XF(1) 

FLCH*PUFF(1) 

XDAs (1 ... 0 1 »FLCH ) »CHORD 
YDAa . 0 1 »PUFF( 2 ) »CHORD 
ARCLs . 0 1 •PUFF( 3 ) •CHORD 
DBFLG=PUFF(4) 

DLTsDLTRfDEFLG 
DLTUsDLTUR'DEFLQ 
DBFLsDEPLG*BOGEN 
ARCLUa . 0 1«PUPF( 5 ) •CHORD 
MSPLIal 

CALL FLAP ( XF , YF , BETAF , NQRS , XDA , YDA , ARCL , ARCLU ,DEFL,X , Y ,XP > NQ ) 
GO TO 176 
C ALFA CARD 

333 IF(NUPA.EQ.O) GO TO 330 
HOMAGaNUPA 
AGAM(2)«FL0AT(NUPE) 

IP(NUPA.EQ.I) AGAM(10)aFL0AT(NUPE-».1) 

330 IP(NUPU.EQ.O) GO TO 335 
DO 331 Ia1,14 

331 ALCA(l)aPUFF(I) 

328 ITIT1aNUPI/2+l 
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ITIT2sNUPI-2»ITIT1+3 POOR QUALITY 

DARFsO. 

IF(ITIT2.NE. 1)DARF=1. 

NALsIABS(NUPU) 

IF(NAL.GT. 14)NAL=4 

335 DO 334 1=1, NAL 
PA=ALCA(I) 

IF(PA.LE.-99.) PA=RS(30+I) 

IF(PA.GT.-99.) PA=PA+DARF«DARG 
334 ALV(I)=PA 

CALL MOMENKX ,Y ,NQ,XDA,YDA,DEFLG ,M0MAG ) 

IF(AGAM(2).EQ.O.)GO TO 11 
NZF=NQ+3 

IF(ITP.EQ.2.AND.AGAMOO).EQ.1.)NZF=0 
CALL DIA(X,Y,NQ,THK) 

THKP=100.«THK 
D0341 N=1,NQ 
ND=N-1 
XDR=X(N) 

YDR=Y(N) 

DO 340M=1,NAL 
V(M)=ABS(VPR(N,M)) 

VQs1.-V(M)»V(M) 

IF(ITIT1.EQ.2)V(M)=VQ 

340 CONTINUE 
NZT=NZPZ(1,NZF) 

IF(NZT.NE.NNESE.AND.N.NE.1)GO TO 341 
NZFsO 

DO 339 Ms 1, NAL 
339 P(M)sALV(M)-DARP»DARG 

332 IF (ITP . EQ . 1 ) WR IT E( IDRU , 337 ) NZT , NUPRO , THKP , ( P (M ) ,Ms 1 , NAL) 

337 FORMAT (A1 ,8HAIRF0IL ,I4,F8.2, 1H>,F9.2, 13F8.2) 
IF(ITP.EQ.2)WRITE(IDRU,336)NZT,NAMP, THKP, FLCH,DEFLG,(P(M),M=1, NAL) 

336 FORMAT (A1 ,8HAIRF0IL , 12A1 ,F8.2, HHJI THICKNESS, PI 0. 2, 6H» FLAP, 
•F8.2,19H DEGREES DEFLECTION/2 3X , 14f8. 2) 

IF(ITP.EQ.2)NZT=NZPZ( 1 ,0) 

NZT=NZPZ(1,0) 

WR ITE( IDRU ,338) NZT ,CPV( ITIT 1 ) , ( ALTX (M , ITIT2 ) ,Ms I , 4 ) 

338 FORMAT (A1,3H N,7X, 1HX,8X, 1HY,5X,A9,60HDISTRIBUTI0NS FOR THE ABOV 

»E ANGLES OF ATTACK REUTIVE TO THE,4A4) 

NZT=NZPZ(1,0) 

341 WRITE(IDRU,342)ND,XDR ,YDR , (V(M) ,Ms I ,NAL) 

342 FORMAT (I4,F10.5,F9.5, 14F8. 3) 

IF( ITP. EQ. 2)00 TO 11 
NZTsNZPZ(1,0) 

WRITE( IDRU , 344) NZT ,DARG ,CM ,ETA 

344 FORMAT (A1,8HALPHA0 s,F5.2,8H DEGREES, 3X,5HCM0 s,F7.4,3X, 

•5HETA s,F6.3) 

GOT0 11 

STOP 

END 
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SUBROUTINE TRAPRO 

DIMENSION PLS(2) ,FLA(2) ,DRAK(2) ,DRAM(2) ,AC(5,4) ,D(4) , 

1WSK2) ,WCIC2) ,FINT(3) ,A(£) ,HK(2) ,R(3) ,FKERN(30) 

DIMENSION CFP(2) ,AKK(2) ,PHIS(2) ,PHIW(2) ,AKN1(2) ,F(2) ,G(2) 

DIMENSION XRT(2),YRT(2) 

C0M10N/ EA/ILBS , IDRU , ISTA , NNESE 

C0MM0N/PRAL/DLT,DLTU,ALN,ALV(14),NAL,ITP,NAMP(12) ,CML(14) ,CRL(14) 

1 ,CPV(2) ,ALTX(4,2) ,DARF,ITIT1 ,ITIT2 

COItWN P1(121),P(121),XP(121),YP(121),PUFF(14),AGAM(14),X(121) , 
1Y(121) ,DS( 122) ,VF( 121) ,ARG(121) ,ANI(28) ,ALFR(29) ,IZ2,KFU,NQ,NUPRO, 
2 JAB , JST , CM , ETA , ABFA , PI , BOGEM , DARG , PURES( 1 3 ) .G AP ( 4 50 ) , ALFA ( 29 ) 

DATA ABSZ/0./ 

CALL WANDEL(NUPR0,NAMP,12,5) 

ALFR(JAB<»-1)sO. 

ABZTsANKJAB) 

IF(ABS(ABZT-ABSZ).LT. .1) GO TO 14 
IBsINT(.25»ABZT+.1) 

MQs2«IB 

NKRs2«MQ 

ABSZsFLOAT(NKR) 

ABGRa360./ABSZ 
HABGRa.5*ABGR 
PURES(8)sO.O 
DO 8 Mai, IB 

ARIaFL0AT(MQ<t>U2»M)«HABGR 
8 FKERN (M ) a ABGR*COSG ( ARI ) / ( SING ( ARI ) *PI ) 

14 MAGAMaINT(AGAM(3)) 

NQaNKRKi 

IF(MAGAM.EQ.O) GO TO 22 
NZTaNZPZ(3,0) 

WRITE(IDRU,82)NZT 

MCTaO 

22 DO 23 Ia1,29 

23 ALFA(I)aALFR(I) 

Ia1 

Jal 

24 FLS(J)a PURBSC I ) *ABFA 

4 CALL DRAW(WC,WS,WL,.6,.1.,FLS(J),ABGR,1) 

CALL DRAU(WCI(J),USI(J),ULI,-.6,.1.,FLS(J) ,ABGR,1) 

WCI(J)a UCI(J)i>WC 
USI(J)a WSI(J)4>WS 
ULI a ULI«WL 
C4aWLI 
CSa.VrLI 

FU(J)a PURES(I-i>1) • ABFA 
IF(PU(J))25,25,26 

25 DRAK(J)a 0 
DRAM(J)a 1. 

GOTO 34 

26 1P(J.EQ.2) 00 TO 401 
UI a COSa(ABQR*PLA(J)) 
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IF(PURES(I+2)-1.)27,30,29 OF POOR QUALITY 

27 DRAK(J)= .1»PURES(I+3) 

28 DRAM(J)= . l»PURES(I+4) 

GOTO 34 

29 DRAK(J)=(( .1*PURES(I+4))»«(-10./PURES(I+3))-1.)*(1-+WI)/(1.-WI) 

DRAK(J)s RUND(DRAK(J),1000.) 

DRAM(J) = . 1»PURES(I+3) 

GO TO 34 

30 AA = .05»(1— WI )*PURES(l+3) 

WILN = AL0G(.l*PURES(I+4)) 

FMIT = .5 

MIT s 0 

31 FM = -WILN/ALOG(AA/FMIT +1.) 

MIT=MIT+1 

IF(ABS(FM-FMIT)-1.E-6) 33,32,32 

32 FMIT = FM 
GO TO 31 

33 DRAM(J) = RUND(FM,1000.) 

DRAS s .05*PURES(I+3)»(WI+1.)/FM 
DRAK(J) s RUND(DRAS,1000.) 

GO TO 34 

401 DRAK(J)s.1»PURES(I-t-3) 

DRAM(J)s1.0 

34 Is 1+5 
Js J+1 

IF(J-3)24,38,38 iS 

38 MER s 0 
WSl(2)s -WSI(2) 

IWPPMsINT(PURESdO)) ~ 

ITM0DsINT(PURES(11)) 

RUFsIOO. 

IF(ITM0D.GE.4.AND.ITM0D.LE.6)RUFs1000. 

ITMRsITMOD 

SHKS s .1»PURES(12) '• 

HKSTs.1»ABS(PURES(13)) 

35 DO 36 Js1,4 : 

36 AC(1,J)s 0. 

ALPA(JAB)sALFA(1) 

ALIV s 0. :• 

SINAI s 0. ; ' 

COSAI s 1. ■ 

FNI s 0. 

TNAIsO. ‘1 

Js1 i! 

37 CSAIP s C0SG(2.»ALFA(J)) 

SNAIP s SINO(2.*ALFA(J)) [’ 

TNAIPsTNG(ALFA(J)) d 

IF(J-JST-1)40,39,40 

39 AC(2,1)s SINAI 

AC(2,2)s .1..COSAI ; 

AC(2,3)» +1. 

AC(3,1)> -SNAIP 


[!® 
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AC(3,2)s 1.+CSAIP 
AC(3,3)=-1. 

AC(4,1)s COSAl-CSAIP 
AC(4,2)s SINAI-SNAIP 
AC(4,3)= 0. 
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AC(5,1)=0. 

AC(5,2)=0. 

AC(5,3)=0. 

AC(2,4)=-2.»THAI 

AC(3,4)=2.»TNAIP 

AC(4,4)=0. 

AC ( 5 , 4 ) s2 . • ( TNAI-TNAIP ) 

ALIS s ALIV 
ALISP = ALFA(J) 

CMTO 41 

PII 3 CSLG(HABGR*PNI>90.,AL1V) 

PIIPs CSLG(HABGR*FNI-90.,ALFA(J)) 

PBsFNI«HABGR*BOGEN 

AC( 1 , 1)s-FIIP»SNAIP-»FII»SINAI+(COSAI-CSAIP)»PB *AC( 1 , 1) 

ACO ,2)*-FII«( 1 .♦COSAI)-fFIIP«( 1 .♦CSAIF ;+(SINAI-SNAIP)»PB^AC( 1 ,2) 
AC(1,3)»-FIIP ♦ FII ♦ ACO,3) 

JCKsJ-1 

JNIsJ.JAB-1 


IF(JCK.EQ.O.OR.JMI.EQ.O) GO TO 402 
FSIsSNLG(HABGR*FNI) 


AC ( 1 , 4 ) sFSI* ( 2 . »TNAI-2 . •TNAIP ) -2 . »TNAI»FII+2 . •TNA1P«F IIP+AC (1,4) 
402 CONTINUE 

41 IF(J-JAB-1 >42,43,45 

42 ALIV sALPA(J) 

SINAIsSNAIP 
COSAIsCSAIP 
TNAIsTNAIP 
FNI s ANI(J) 

J*J+1 


GO TO 37 
43 Jsl 

IF(FLA(J)) 47,47,49 

49 CALL DRAH(UC,WS,WL,DRAK(J),ORAM(J),FU(J),ABGR,0) 
AC(1,1)s VC-«> AC(1,1) 

AC(1,2) s US 4>AC(1,2) 

AC(1,3)>ULi>AC(1,3) 

47 J»2 
UHCSsO.O 
UUS6«0.0 
C6s0«0 

IP(FLA(J)) 410,410,411 

411 CAa DRAH(1IC,tlS,HL,ORAK(J),1.0,PU(J),ABQR,0) 
WC6«UC 
UVS6«.US 
C6«.WL 
410 OONTINUB 

PHIS( 1 } «PL8( 1 ) •ABOR*BOGBN 
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PHIW( 1 ) =FLA( 1 ) •ABGR»BOGEN 
PHIS(2)3(3dO.-FLS(2)»ABGR)»BOGEN 
PHIW(2)=(360.-FLA(2)«ABGR)*B0GEN 
DO 425 J=1,2 

G(J)sDRAK(J)/(1.0+COS(PHIW(J))) 

F(J)=1.0-G(J)»COS(PHIW(J)) 

IF(G(J)»»2-F(J)»»2) 420,421,422 

420 TABG=((F(J)-G(J))«TAN(PHIW(J)/2.) )/SQRT(F(J)»»2-G(J)»»2) 

TARG 1 sATAN ( TARG ) 

CFP(J)=-((4.»(G(J)-F(J)))/SQRT(F(J)«»2-G(J)«»2))»TARG1 
GO TO 425 

421 CFP(J)=0.0 
GO TO 425 

422 GFLN2=SQRT(G(J)«»2-F(J)«»2) 

GFLN1=(G(J)-F(J) )»TAN(PHIW(J)/2.) 

GFLN=ALOG(ABS( (GFLN1-KJFLN2)/ (GFLNUGFLN2) ) ) 

C0EFF=-(2.«(G(J)-F(J) )/SQRT(G(J)«»2-F(J)*»2) ) 

CFP(J)aCOEFF»GFLN 

425 CONTINUE 
SRMa.6 
AK=1.0-SRM 
BKsl.O-^SRM 
DO 450 Js1,2 

AKN 1 ( J ) s SRM/ ( 1 . 0-COS ( PH IS ( J ) ) ) 

IF(BK-2.»AKN1(J)) 430,431,432 

430 AKSRsSQRT(2.«BK»AKN1(J)-BK»»2) 

AKT1s2.*AKN1(J)-BK 

C0EFKsAKT1/(2.»AKSR) 

AKLN 1 sAKT 1 »TAN ( PHIS( J ) / 2 . ) 

AKLNs ALOG ( ABS ( ( AKLN 1 -t-AESR ) / ( AKLN 1 - AKSR ) ) ) 

AKK(J)aCOEFK*AKLN 
GO TO 450 

431 AKK(J)sO. 

GO TO 450 

432 AKSR1=SQRT(BK*»2-2.«AKN1(J)«BK) 

AKT2a2.»AKNl(J)-BK 

AKT3aBK>2.«AKN1(J) 

AKT AN a ( AKT 3*TAN ( PH IS ( J ) / 2 . ) ) / AKSR 1 
AKATANaATAN(AKTAN) 

AKK ( J ) a ( AKT2/ AKSR 1 ) *AKATAN 
450 CONTINUE 

AC1a-2.»DRAM(1)«PHIW(1) 

AC2aDRAM(l)«CFP(1) 

AC(1,4)>AC(l,4)fACUAC2 
C 1 AaSQRT ( ( AK4-2 . *AKN U 1 ) ) / AK ) 

C1Ca4.»C1A»ATAN(C1A»TAN(PHISO)/2.)) 

C1«4.0«PHIS(1)-C1C4.4.0»AKK(1) 

C2AaSQRT( (AK4.2.»AXN1 (2) )/AK) 

C 2C « 4 . »C2A*ATAN ( C2A*TAN ( PHIS( 2 ) / 2 . ) ) 
C2a8.0»PI-4.0«PHIS(2)>C2C-4.0»AKK(2) 

C3«2.0»PHIH(2)-CPP(2)-4.0»PI 

D{l)aH5I(1)«(C3»C5-C2«C6)-HSl(2)«(C3*C4-C1»C6)+i«S6»(C2*C4-C1«C5) 
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D(2)S-(WC1(1)»(C3«C5-C2*C6)-WCI(2)»(C3»C4-C1*C6)+WWC6*(C2*C4-C1«C5 

D(3)--WCI(1)»(C3»WSI(2)-C2»WWS6)-WCI(2)«(C3»WSI(1)-C1»WHS6)fWWC6»(C 

12»WSI(1)-C1«WS1(2)) 

D(4)s-(WCin)«(c6»WSI(2)-C5»WWS6)-WCI(2)»(C6»WSI(n-C4«WWS6)+WWC6» 

1(C5»WSI(1)-C4»WSIC2))) 

A(1)sD(1)«AC(1,1)+D(2)»AC(1,2)^D(3)»AC(1,3)+D(4)»AC(1,4) 

A(2)xD(1)»AC(2,1)+D(2)»AC(2,2)+D(3)»AC(2,3)fD(4)»AC(2,4) 

A(3)=D(1)»AC(3,1)+D(2)«AC(3»2)+D(3)*AC(3,3)+D(4)»AC(3,4) 

A(4)sD(1)«AC(4,1)+D(2)«AC(4,2) 

A(5)=D(4)»AC(5,4) 

C SOLUTION OF TRANSCENDENTAL EQUATION 
53 IsO 

FV » 9.E9 

PHISH s .5 •(ALIS^ALISP) 

60 CSLI 3 CSLG( PHISH, ALIS) 

CSLIPx CSLG( PHISH, ALISP) 

SNLIxSNLG(PHISH-»-90. ) 

FP xA ( 1 ) ♦A ( 2 ) •CSLI+A ( 3 ) •CSLIP+A ( 4 ) •BOGEN* ( 90 . ♦PHISH ) ♦A( 5 ) *SNLI 
IF(I.GE.20) GO TO 66 
IF(ABS(FP)-ABS(FV).LT.-.5E-9) GO TO 62 

1x20 

PHISH 3 PHISH - PDIF 
GO TO 60 

62 PDIF3-FP/(A(2)/(PHISH-AL1S)>(A(3)/(PHISH-ALISP))^A(5)«PHISH) 

I 3 1+1 

65 FVxFP 

PHISH 3 PHISH ♦ PDIF 

IF(PHISH.LT.ALIS. AMD. PHISH.GT. ALISP) GO TO 60 
WR nE( IDRU , 64 )mER , ITHOD 

64 FORMAT (66H0TRANSCBNDBNTAL EQUATION HAS DIVERGED. CHECK TRA1 AND 
•TRA2 CARDS., 12H CTERATION , 12 , 8H MODE ,11) 

STOP 

66 ANI(JST) 3 (PHISH4-90.)/HABQR 

A JCx AC ( 1 , 1 ) kAC ( 2 , 1) •CSLIfAC (3,1) *CSLIP^AC (4,1) •BOOEN* ( PHISH-»-90 . ) 
AJSmAC( 1 ,2)-»>AC(2,2)*CSLI^AC(3,2)*CSLIP^AC(4,2)*BOGEN*(PHISH4>90. ) 
AJTxAC( 1 ,3)'*'AC(2,3)«CSLI+AC(3,3)»CSLIP 
AJ43AC(1,4)^AC(2,4)«CSLl4.AC(3,4)*CSLIP^AC(5,4)*SNLI 
DDxNCI ( 1 ) • (NSI ( 2 ) »C6-WS6«C5 )-NCl ( 2 ) • (NSI ( 1 ) •C6-WS6»C4 ) ♦iftIC 6» ( WSI 
1(1)»C5-NSI(2)»C4) 

DRAM( 2 ) * (UCI ( 1 ) • ( A JS«C5-A JTmSI ( 2 ) ).UCI ( 2 ) • ( A JS«C4-A JT«WSX ( 1 ) ). A JC 
1 • (NSI ( 1 ) •C5-NS1 ( 2 ) *04 ) )/DD 

69 HK(1)>(.AJC*afSI(2)«c6.NHS6«C5)-WCI(2)«(AJT^S6.AJS*C6)^tiWC6«(AJT 
1*H81(2).AJS«CS))/DD 

HK(2)i(NCl(1)*(AJT«WS6.AJS«C6)^AJC*(H8I(1)*C6.tNS6*C4)^WC6«(AJS* 
1CMI8I(1)*AJT) )/DD 
HK8 « NK(1)<HK(2) 

C 

HPPWD«( 1 .OfDIUK( 1 )♦( 1 .0.D1UX( 1 ) )*008(PHIW( 1 ) ) ) 

HFPUU*(DIUM( 1 )*ORAX( 1 ) )/VPPW)D 
NPPSUsd . 125*NK( 1) )/ ( 1 .0.^8(PHZ8( 1) ) ) 

UPPUsNPPUU^WPPau 
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MPPWLDsd .0*DRAK(2)+( 1 .0-DRAK(2) )«COS(PHIW (2) ) ) 

WP r ;l = ( dr am ( 2 ) »orak ( 2 ) ) / wppwld 

WPPSL=( 1 . 125»HK(2) )/(1 .0-C0S(PHIS(2) ) ) 

WPPL=WPPWL+WPPSL 

IFdTMOD.EQ.O.OR.ABS(HKS-SHKS).LT.HKST) GO TO 74 
IF(MAGAM.LT.2.AND.HAGAM-MER.NE. 1) GO TO 100 
GO TO 76 
74 ITMOD =0 

IF(MAGAM.EQ.O) GO TO 300 

76 NZTsNZPZ(2,JAB+4) 

WR ITE( IDRU , 77 ) NZT , NUPRO ,MER , ITMR 

77 FORMAT (A1 ,42HTRANSCENDENTAL EQUATION RESULTS AIRFOIL ,14, 

•12H ITERATION, 12, 8H MODE ,11) 

NZTsNZPZ(1,0) 

WRITE(IDRU,78)NZT 

78 FORMAT (A1,69H NU ALPHA* OMEGA' OMEGA K MU KH 

• LAMBDA LAMBDA*) 

JHs 1 

DO 85 JN=1,JAB 

79 NZT*NZPZ(1,0) 

IF(JN.NE. 1.AND.JN.NE.JAB) GO TO 83 
XI = .5*(1.-*' COSG(FLA(JH)*ABGR)) 

WHK = (1.>DRAK(JH)*(1.-X1)/X1)**(-DRAM(JH)) 

USTRs DRAM(JH)*DRAK(JH)/X1 

WRITE(IDRU,82)NZT,ANI( JN) ,ALFA( JN) ,WSTR,WHK,DRAK( JH) ,DRAH(JH) ,MX( J 
1H),FLA(JH),FLS(JH) 

82 FORMAT (A1,F6.2,F8.2,F8.3,F7.3,F9.4,F9.4,P8.3,F8.2,F7.2) 

JHs2 

GO TO 85 

83 WRITE(IDRU,82)NZT,ANI( JN) ,ALFA(JN) 

85 CONTINUE 

WR ITB( IDRU , 2000) WPPU ,H?PL, IHPPM ,MCT 
2000 FORHAT(/5X,5HWPPUs,F7.3t5X,5HUPPLi,F7.3.5X,11HWPP ITH00E«,I2,3X, 
210HITERATIONs,I2/) 

IF(ITMOD.EQ.O) GO TO 300 
100 IF(MER)103,102,103 

102 DAL > .1 
CO TO 104 

103 IF(HKS-HXSV.EQ.O.)GO TO 74 

DAL > (SHKS.HKS)*DAL/(HKS.HKSV) 

DALDsDAL 

IF( ITP.EQ.O)OAL>RUND(DAL,RUF) 

IF(MACAH.EQ.O.)GO TO 1004 
NZTaNZPZ(2,0) 

WRnEClDRU, 1003)NZT,HER,HKS,DAL0,DAL 

1003 FORMAT (A1,10H ITERATION, 12, 3X, SHE S ■,F9.6,3X,7HOELTA «,F12.8, 
*3X,9HROUNOED *,P6.2) 

1004 IF(DAL.EQ.O.)GO TO 74 
1F(MER.0E.3*AND.ABS(DALV).LE.ABS(DAL))00 TO 74 
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104 DALVsDAL 

IFaTMOD.OE.4)GO TO 113 
DO 111 Jsl.JAB 

1P( ITMOD . ME . 2 . AND . J . LE . JST ) ALFA ( J ) a ALFA ( J ) 4-DAL 
IFdTMOD.HE.I.AAlD.J.GT.JST) ALFA( J)sALFA( J)-0AL 


111 CONTINUE 
GO TO 112 

113 IP( ITMOD .OE. 7)00 TO 114 

IF( ITMOD . NE . 5 )DRAK ( 1 ) *DRAK ( 1 ) ♦DAL 
IFv ITMOD . NB . 4 )DRAK ( 2 ) sDRAK ( 2 ) kDAL 
GO TO 112 

114 IP(ITM0D.NE.8)ALFA(JST)^ALFA(JST)fDAL 
IF(ITM0D.NE.7)ALFA(JST>1)»ALFA(JST*1)-DAL 

112 HKSVsHKS 
MER aMER>1 
GOTO 35 

C 

300 OHPPaWPPO-WPPL 
CFUNCsABS(DUPP) 

NPPTOLa.001 

IF(CPUNC.LE.WPPTOL) GO TO 301 
IF(IWPPM.EQ.O) GO TO 301 
IF(MCT.GT. 15) GO TO 301 
IF(MCT.NE.O) GO TO 19 
DBU.1 
GO TO 20 

19 DBU-DHPP*DBL/(DMPP-DWPPV) 

IF(DBL.EQ.O.) GO TO 301 

20 IF(IHPFM.BQ.4) GO TO 7 
IF(IHPPM.EQ.3) GO TO 6 
IF(IHPPM.EQ.2) GO TO 9 

15 DRAK(1)sDRAK(1)t>DBL 
GO TO 21 

9 DRAK(2)«DRAK(2)oDBL 
00 TO 21 

6 DRAH(1)aDRAM(1)i>DBL 
00 TO 21 

7 FLSH<FLS(2)-«>DBL 

IF (FLSN.LT.0.0.0R.FLSN.0T.FU(2)) GO TO 12 

FLS(2)«FLSN 

DBLVsDBL 

MfPPV-OWPP 

HCT«MCT4>1 

I«6 

00 TO 4 
12 IHPPMbI 
DCU.1 
00 TO 15 

21 OBLV«DBL 
MfPPV«OIIPP 
HBRiNBII^I 
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IF(ITMOD.EQ.O) MER=0 

MCT=MCT+1 

GO TO 35 

301 AX1 = .5»(C0SG(PHISH-ALFA(JST+1))/SING(PHISH-ALFA(JST+D) 

1 -COSG(PHISH-ALFA(JST))/SING(PHISH-ALFA(JST))) 

AKP =AK1*180./9. 8696044 
PHIM = 0. 

NU=1 
1 = 1 
ANU =0. 

JH=0 
VI= 0. 

302 JH=JH+1 

FF1 = 'X)SG(ABGR»FLA(JH)) 

FF2 = jRAK(JH)/(l.+FFl) 

FG1 = COSG(ABGR»FLS(JH)) 

FG3 = .6/(FG1-1.) 

304 Vis VI - CSLG(PHIM-90,, ALFA(I)) 

GO TO 310 
306 ARGN = ANU 

IF(ANU.GT..5* ABSZ)ARGNs ABSZ - ANU 
CSP = COSG(ARGN»ABGR) 

FIsO.O 

IF( ARGN . LT. FLA( JH ) )F1 sDRAM( JH ) »ALOG ( (CSP-FF1 ) •FF2+ i . ) 

GIsO. 

IF(ARGN.LT.FLS(JH))G1=-HK(JH)«AL0G(1.-((CSP-FG1)«FG3)»*2) , i 

P(NU)S F1+G1+CSLG(ANU»HABGR-90.,ALFA(D) ♦ VI 
P1(NU)sP(NU)-AK1«ABS(SING({ANU«HABGR - 90.) - PHISH)) 

NU 3 NU ♦ I 
ANUs ANU-»- 1. 

310 IF(ANU-ANKI) >306,306,312 
312 IF(ANU- ABSZ) 3 14, 320, 320 
314 PHIM s ANI(I)»HABGR 

VIsVI+CSLG(PHIM-90. ,ALFA(I) ) 

I s 1+1 • ; 

IF(I-1-JST)304,302,304 : ; 

320 PSsO. 
e2s0. 

DO 324 Is1,NKR 
PSsPS+P(I) 

BI s 2«(I-1) 

324 B2 s B2 ♦ SING(BI»ABGR)»P(I) ' 

VI s 2.»EXP(PS/ABSZ) 

SXI 3 .00000000 

SYsO. - ; 

D0328 Ns1,NQ 

QsO. 

D0326 Ms1,IB 

MN s N ♦ 1 + MQ - 2»M Tl 

MM s 2«N . MN 

IF(MN.GT.NKR) MN s MN - NKR 

• ] 

i i 
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C 


IF(MM.LT.I) MM = MM ♦ NKR 
326 Q = Q+ FKERN(M)*(P1(MN)-P1(MM)) 
ANU= N-1 


ZP = ANU*HABGR - 90. 

ZL = COSGCZP - PHISH) 

ZL s ABS((1.-ZL)/(1.+ZD) 

IF(ZL.NE.O.)ZL3ALOG(ZL) 

ARG(N) = Q - AKP«SIMG(ZP-PHISH)»ZL ZP 
VF(N) = V1»EXP(-P(N)) 

WV = COSG(ZP)/VF(N) 

XP(N)= WV»SING(ARG(N)) 

YP(N)s-WV«COSG(ARG(N) ) 

SXIs SXI+ XP(N) 

328 SY s JY ♦ YP(N) 

SX : SXI 

XPK s SX/CABSZ >1.) 

YPK s SY/(/.BSZ -1.) 

DO 329 N:2,NKR 
XP(N}:XP(N)-XPK 

329 YP(N)=YP(N)-YPK 

CALL CINT(XP,X,NQ,IZZ) 

CALL CINT(YP,Y,NQ,IZZ) 

RQV = 0. 

DO 330 Ns2,NKR 
RQaX(N)»X(N)+Y(N)»Y(N) 

IF(RQ.GT.RQV)LsN 

330 RQV = RQ 

DO 327 I a 1,3 
lEPPL a L-2+I 

327 R(I)aSQRT(X(IEPPL)«X(IEPPL)-*>Y(IEPPL)»Y(IEPPL) ) 

333 TAU a (R(3)-R(1))/(4.»(R(2)+R(2)-R(1)-R(3))) 

XNAS a X(L)fTAU»(X(L+l)-X(L-1)+2.»TAU»(X(L+1)+X(L-1)-X(L)-X(L))) 

YNAS a Y(L)+TAU«(Y(Lfl)-Y(L-1)+2.»TAU«(Y(L>1)<*-Y(L-1)-Y(L)-Y(L))) 

SO a XNAS*XNAS <•> YNAS«YNAS 

ATaXNAS/SQ 

Ba YNAS/SQ 

STREP a 1./SQRT(SQ) 

ETA a ABSZ*STREF/PI 
CM a .5»ETA»STREF»B2 

DARG a 19.09859 •(3.*YNAS/XMAS - (YNAS/XNAS)*»3) 
1F(ABS(SX)4>ABS(SY).LT..0001«ABSZ) GO TO 335 
SXaSTREF«SX»200. 

SYaSTREF«SY«200. 

NZTaNZPZ(2,0) 

WRHE (IDRU,334) NZT,SX,SY 

334 FORMAT (A1,14HWARNING - SX a,F6.3»3X,4HSY a,F6.3) 

335 CONTINUE 


IF(IWPPM.EQ.O) GO TO 605 
CHORDa ( 4 . 0'-'PI ) / ( ABSZ»STREF) 
PHIVr 1 aCOSQ ( FU( 1 ) •ABOR ) 
PHIS 1 aCOSO ( PLS ( 1 ) *ABGR ) 
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PHIWF= 1 . 0+DRAK ( 1 ) »( (1 . 0-PHIW 1 ) / ( 1 . 0+PHIW 1 ) ) 
W0=(PHIWF»»(-DRAM(1)))»(.64)»»HK(1) 

APOs-ALOG(V1 )-ALOG(WO)+ALOG(2.0) 
a0=(ARG(1)+90.)«(PI/l80.) 

ALFA1=ALFA(1)»(PI/180.) 

DELTE=.5»(ATAN(Y(2)/X(2))+ATAN(Y(NKR)/X(NKR))) 

GAMMA:ALFA1-DELTE 

605 DO 331 Ns2,NQ 
XR=X(N) 

X(N)= 1-B«Y(N)-AT»XR 
Y(N)= B»XR -AT*Y(N) 

ARG(N) = ARG(N) - DARG 

WQ = (XP(N)+XP(N-1)-XPK-XPK)»»2 + (YP(N)+YP(N- 1 )-YPK-YPK)»»2 

331 DS(N-I) s STREF»SQRT(WQ)«(1.+.6666667 »((XP(N)»yP(N-1) 
1-XP(N-1)»YP(N))/WQ)*»2) 

NHKW:NQ/12 

DLT = Y(NHXW)/(BOGEN*(1.-X(NHKW))) 

NHKU:NQ-NHKW-*-1 

DLTUs-Y ( MHKW ) / I BOGEN* ( 1 .-X ( NHKW ) ) ) 

X(1) = 1. 

332 ARG(l) = ARG(1) - DARG 
346 ITPsI 

IFdWPPM^EQ.O) GO TO 606 

XTsXd) •CHORD 

YT=Y(1) •CHORD 

XRT(1)=X(2)^CH0HD 

XRT(2)=X(NKR)^CH0RD 

YRT(1)sY(2)^CH0RD 

YRT(2)sY(NKR)^CHORD 

CALL RCAL(XT,YT,XRT,YRT,AP0,Q0,ALFA1 , PXT , PYT , NKR ) 
PST=PXT^COS(GAMMA)+PYT^SIN(GAMMA) 

PNTs-PXT^SIN (GA^f^A) +PYT»C0S(GAMMA ) 

CPST=2.0^CHORD^PST 

CPNTs2.0^CH0RD^PNT 

WRITE(IDRU,602) 

602 F0RNAT(”0",11X,45HTHE TRAILING EDGE PRESSURE GRADIENT IS FINITE) 
WRITE( IDRO ,603) ALFA( 1 ) 

603 FORMAT(15X,25HWHEN THE ANGLE OF ATTACK=,F4. 1 , 1X,7HDEGREES) 
WRnE(IDRU,604) CPST.CPMT 

604 F 0 RMAT( 8 X, 19 HIN THAT CASE, CPST»,F10.3,7X,5HCPNTa,F10.3) 

606 ALNsDARG 

IF(P0RBS(13).GE.0.)G0 TO 11 
PURES(12)a10.^HKS 
PURES(13)s>00001 
11 RETURN 
END 
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FUNCTION TNG (A) 

TNG=SING(A)/COSG(A) 

RETURN 

END 


FUNCTION ^LG(A) 
SNLG=ALOG(ABS(SING(A) ) ) 
RETURN 
END 
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SUBROUTINE RCAL(XT,YT,X,Y,AP0,Q0,ALFA1 , PXT , PYT , NKR ) 
DIMENSION PX(4) ,PY(4) ,X(2) ,Y(2) ,THETA1(2) 

COMPLEX CEXP , CMPLX ,CON JG , ZETA , EIA , EMIA , ZTEXP , 

1 ZETA2 , WZ , ZT2P , ZT3P , WTP , WT2P , ZETA 1 , Z , ZT4P , ZT , 

2ZT2PC , PTCAL , C ( 2 ) , ZETA 1 P( 2 ) , WT3P , WTPC , PTCOEF 
PI=3. 141592653589 
ANKRrFLOAT(NKR) 

THETDEG=360./ANKR 
THETAK 1)=THETDEG 
THETA 1 ( 2 ) - j 60. -THETDEG 
ZT2XP=CMPLX(AP0,Q0) 

ZT2P=CEXP( ZTEXP) 

COSAsCOS(ALFAI) 

.JINA=SIN(ALFA1) 

3IA=CMPLX(C0SA,SINA) 

SMIA=CMPLX(COSA,-SINA) 

WTP=(EIA+EMIA) 

WT2P=(-4.0»EIA-2.0«EMIA) 

WTPCsCONJG(WTP) 

WT3P =6 . 0» ( 3 . 0»EIA+EMI A ) 

ZT3P=ZT2P»(WT2P/WTP) 

ZT=CMPLX(XT,YT) 

ZT4P=0.0 

M=1 

DO 20 1=1,2 
Z=CMPLX(X(I),Y(I)) 

THETAs THETA Kl)»( PI/180.) 

XI=COS( THETA) 

ETA=SIN( THETA) 

ZETA=CMPLX(XI,ETA) 

ZETA1P(M)s(ZETA-1.0) 

C (M ) = ( 1 20. /ZETA 1 P (M ) »»4 ) • ( Z-ZT- . 5»ZT2P«ZETA 1 P(M ) «»2 
1-(ZT3P/6.)»ZETA1P(M)*»3) 

M=2 

20 CONTINUE 

ZT4P=. 10»(C( 1 )+C(2) )+ZT4P 
ZT2PC=C0NJG(ZT2P) 

PTCOEFs- (WTPC )/ ( 3 • 0»ZT2PC«ZT2P»»3) 

PTCALsCON JG ( PTCOEF* ( ZT2P*WT3P-ZT4P*WTP ) ) 

PXT = REAL (PTCAL) 

PYTsAIMAG(PTCAL) 

RETURN 

END 
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igure 1. Form of the velocity distribution for mexlmum lift on a 
single-element airfoil. 
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Figure 3. Kennedy and Marsden airfoil potential flow velocity distribution 
(Eppler panel method code) compared with viscous analysis 
(GRUMFOIL code) and experimental results [18]. ALPHA ■ 4.2 DEG 
(relative to chord-line), RE • 1 x 10^. ,( 
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Figui'e 7. Example von Mises airfoil generated using a 
four-term transformation. 
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Figure 9. Example von Mises airfoil generated using a 

six-term transformation and a mapping circle radius 
equal to the minimum allowable for a physically 
realizable airfoil. 
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Figure 10. Symmetrical von Mises airfoil having real transformation 
coefficients equal to those of the airfoil shown in 
Figure 7. 
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Figure 11. Von Mises airfoil generated with a mapping circle 

radius less than the minimum allowable for a physically 
realizable airfoil. 
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Figure 12. Von Mises airfoil generated with a mapping circle 
radius greater than the minimum allowable for a 
physically realizable airfoil. 
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Figure 13. Flat plate airfoil generated by simplifying the 
von Mises transformation to that of Joukowsky. 
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Figure 14. Von Mises airfoil resulting from small perturbation 
of zero locations giving the flat plate result. 
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Figure 15. Von Mises airfoil resulting from small perturbation 
of zero locations giving the flat plate result. 
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Figure 16. Von Hises airfoil having largest VTE attainable 
for the given real parts of the transformation 
coefficients. 
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Figure 17. Von Hises airfoil resulting from the small perturbation 
of the zero locations used in generating the airfoil 
of Figure 16. 
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Approxiaat* r«l«tlonship between the BexiMue trailing edge 
velocity ratio and the thickness ratio for reasonable 
syMetrical airfoils. 
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Figure 21. Transformation of a unit circle at angle of attack 
into an airfoil. 
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Figure 22. Von Mises airfoil having finite trailing edge 
pressure gradients, CPXT » - 1.37, CPYT ■ 0.0 
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2-plane 




Figure 23. Segmentation of the interval (0, 2ir) in the circle 
and airfoil planes. 
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Figure 24. Typical pressure recovery distribution. 
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Figure 25. Airfoil obtained using Eppler and Somers code [25] in which 
no iteration is performed for achieving a desired tralllnK 
edge closure angle. ALPHA - 8.0 DEG (relative to aero- 
lift line). 
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igure 26. Airfoil obtained using the modified Eppler and Somers code 
in which the integral constraint required for finite 
trailing edge pressure gradients is satisfied. 

ALPHA ■ 8.0 DEG (relative to sero-lift line). 




Figure 27. 


Airfoil resulting from Eppler and Somers code In which K 
is iterated to achieve a desired trailing edge angle. 
ALPHA *8.0 DEG (relative to zero-lift line). 
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Figure 28. Airfoil having finite trailing edge preaaure gradients 
obtained using nodlfled Eppler and Sonars code. 

ALPHA - 8.0 DEG (relative to sero-lift line). 

CPST • - 7,2 CPNT - - 57.0 
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Pltut* 32. Potential flow velocity diatrlbutlon of Strand airfoil » 

Reference [34], calculated uaing the panel nethod analyela 
capability of the Eppler and Sonera code, ahown along with 
the finite trailing edge preaaure gradient airfoil and ita 
daaign velocity diatribution. ALPHA - 12.0 DEC (relative 
to sero-lift line). 
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Figure 37. 


Velocity distribution at the design point and finite trailing 
edge pressure gradients airfoil based on FX 67-K-150 of 
Wortmann. 

ALPHA » 10.0 DEG (relative to zero-lift line). 
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Figure 41. Local boundary layer form factor, ^^2* equivalent surface source velocity 
distribution to account for displacement thickness. Wortmann FX 67-K-150. 
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distribution to account for displacement thickness. Airfoil based 
FX 67-K*150 but having finite trailing edge pressure gradients. 
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Figure 45. Off-deeign pressure dlstrlbuclon from GRUMFOIL analysis. Airfoil based on Wortmann 
FX 67-K-lSO but having finite trailing edge pressure gradients. CL “ .7019 
CD - .0055 CM - -.1211 CPTE = .2019 ALPHA - 0.98 RE - 2x10^ 
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